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1)  It  measures  In  "real  time”,  directly  a  component  of  rotation  rate,  2)  It  Is 
a  local  sensor  of  absolute  rotation,  3)  It  can  basically  be  adjusted  to  be 
most  sensitive  for  slow  events  (daily-monthly  changes  of  the  earth  rotation 
vector)  or  fast  events  (microseisms,  earth  quakes,  local  disturbances  of  earth 
surface),  4)  since  It  Is  a  local  method,  the  beam  can  naturally  be  circulated 
In  vacuum,  atmospheric  disturbances  are  eliminated  from  the  beam  path. 
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Foreword 


The  project  of  constructing  a  large  laser  ring  is  an  end  result  of 
several  questions,  the  more  important  of  them  being:  a)  Why  should  one 
do  it?  (motivation,  expected  benefits,  comparison  with  competitive 
methods),  b)  Is  the  project  feasible?  (sensitivity,  costs, 
compatibility  of  expected  results  with  a))  c)  What  is  the  design? 
(technical  parameters,  stability  considerations,  output). 

Ample  literature  is  available  on  a).  A  recent  survey,  although 
incomplete,  is  ref.  1.  It  had  become  clear  that  the  project  would 
become  worthwhile  if  rotation  sensitivities  of  better  than  lCT^n^ 
(fyr=‘'earth  rate"  =  2ir/day)  can  be  reached.  In  this  case  it  can 
potentially  surpass  in  sensitivity  not  only  the  most  advanced  techniques 
for  measuring  absolute  rotation  of  the  earth,  i.e.  Lunar  ranging, 
Lageos,  VLBI  (ref.  2),  but  it  offers  exciting  possibilities:  1)  it 
measures  in  "real  time",  directly  a  component  of  rotation  rate,  2)  it  is 
a  local  sensor  of  absolute  rotation,  3)  it  can  basically  be  adjusted  to 
be  most  sensitive  for  slow  events  (daily-monthly  changes  of  the  earth 
rotation  vector)  or  fast  events  (microseisms,  earth  quakes,  local 
disturbances  of  earth  surface),  4)  since  it  is  a  local  method,  the  beam 
can  naturally  be  circulated  in  vacuum;  atmospheric  disturbances  are 
eliminated  from  the  beam  path. 

The  extrapolation  to  an  expected  sensitivity  of  <10"^a^  is  done 
from  data  presently  available,  of  the  order  of  lO"^  (with  averaging 
times  of  1  day),  with  the  theory  of  laser  noise  which  is  reasonably  well 
in  hand  at  this  time.  Much  of  the  data  on  noise  were  measured  by  the 
author  during  several  summer  joos  in  industry.  The  1/f-noise  becoming 
overseeable  (ref.  3),  we  may  now  be  able  to  predict  dCT^nc  for 


averaging  times  much  less  than  Id.  Technological  improvements  in  making 
high-quality  mirrors  do  continually  push  back  the  limits  due  to  white 


noise  and  1/f  noise.  The  finished  ring  should  become  a  formidable 
sensor  of  £2g<fc),  for  Fourier  frequencies  larger  than  about  a  reciprocal 
year. 

This  report  provides  a  variety  of  technical  details  on  the  basic 
noise  limitations  in  a  laser  ring  and  on  how  to  remove  obstacles  to 
reach  the  required  sensitivities;  some  of  the  basic  calculations  on 
rings  are  collected  as  well  as  several  novel  insights  into  the  setting 
up  are  reported.  The  major  results  have  already  been,  or  are  in  the 
process  of  being,  published. 

As  a  summary  of  the  first  year  of  involvement  with  the  Seiler  ring 


project,  we  state  that  none  of  the  results  have  yet  negated  the  proposed 
performance  of  it. 
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Report 

a 

Summary:  The  results  obtained  can  be  divided  into  three 

categories:  a)  A  feasibility  study  with  focus  on  quantum  noise  and  low 
frequency  noise  b)  Basic  modeling  of  the  ring  with  Gaussian  beam  and  ray 
matrices  c)  Technical  design:  Effect  of  residual  gas  in  ring  on  quality 
factor  and  light  drag,  scanning  of  beam,  effect  of  misalignment  and 
mismatch  of  source  to  ring,  calibration  procedures. 

The  results  show  no  obstacle  yet  to  the  goal  of  achieving  a 
sensitivity  of  rotation  rate  of  better  than  10'9  (earth  rate)  in  rings 
of  60  m2  size.  Such  a  sensitivity  which  corresponds  to  a  change  of 
earth  surface  velocity  of  smaller  than  4  cm/day,  should  surpass  Lunar 
ranging  methods,  Lageos  methods  as  well  as  VLBI-methods  in  accuracy, 
besides  being  a  Veal-time*  observational  method  for  earth  rotation. 


Content 


I.  Introduction 

II.  Contributions  to  basic  theory  of  ringlaser 

II.  A.  Noise  in  ringlasers:  Quantum  noise  and  1/f  noise 

II.  B.  Astigmatic  Gaussian  beam  in  plane  ring:  Stabil ity, curvature 
radius  of  beam,  location  and  size  of  waist(s) 

II.  C.  Ray  matrix  approach  to  square,  plane  ring  with  equal  mirror 
radii:  Stability,  waist  sizes,  spot  sizes  on  mirrors,  design 
for  minimum  spot  size  on  mirrors  (minimum  diffraction,  smallest 
mirror  size),  design  for  circular  spot  size  on  mirrors  (near- 
confocal  cavity  for  best  mode-matching) 

II.  D.  Injection  errors:  Response  of  ring  to  offset  and  tilt  of 
injecting  ray 

II.  E.  Frequencies  of  Hermi te-Gaussian  modes  in  a  ringlaser 

II.  F.  Injection  errors:  Misalignment  and  mismatching  of  an  injected 
Gaussian  beam.  Effects  of  offset,  tilt,  mismatch  in  creating 
Hermi te-Gaussian  eigenmodes  (H.-G.'s).  Power  loss  and  pulling 
of  (fundamental)  Gaussian 

III.  Contributions  to  design  of  ringlaser 

III.  A.  Scanning  of  an  astigmatic  beam,  and  beam  evolution:  New  methods 
to  find  beam  and  evolution  parameters 

III.  B.  Vacuum  in  ring:  Residual  Fresnel  drag 

III.  C.  Vacuum  in  ring:  Effect  on  quality  factor 

III.  D.  Proposed  calibration  methods  for  large  ring:  Light  drag 

through  gas  flow  from  switchable  calibrated  leaks,  tilting 
(rocking)  the  ring  base 


I.  Introduction 
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The  content  of  this  report  is  mainly  intended  to  facilitate  the 
design  of  large  laser  rings  which  act  as  gyros:  They  sense,  via  Sagnac 
effect,  the  absolute  rotation  of  their  structure.  Versions  with  sizes 
up  to  square  meters  are  by  now  well-known  in  industry  as  R(ing)  l(aser) 
G(yro)‘s.  They  have  been  pushed  to  better  performance  ever  since  the 
basic  proposal  of  Rosenthal  ^  was  published.  To  guide  the  reader:  two 
basic  equations  govern  signal  and  noise  of  such  a  gyro: 

Signal:  Af=(4/XL)A*fi  1) 

This  equation  translates  an  absolute  rotation  rate  t )  of  the  ring  into 
a  beat  frequency  Af  (t); 

Noise:  SAf=hfo/Q2P  2) 

which  represents  white  quantum  noise  fluctuations  per  mode.  Hereby,  SAf 

S 

is  the  power  spectral  density  of  the  output  frequency  noise  Af  . 

This  report  can  be  visualized  as  a  study  to  realize  the  promise  of 
equations  1  and  2.  Most  of  the  results  were  obtained  with  a  view  on 
designing  the  Seiler  ring  at  the  AF  Academy,  but  they  are  of  course 
generally  applicable  to  large  rings.  In  certain  cases  calculations  were 
done  on  an  active  ringlaser.  The  equivalence  in  performance  of  active 
and  passive  rings  6  justifies  such  a  procedure,  although  the  design 
calls  for  a  "passive"  ring  at  this  time. 

Whenever  a  certain  topic  treated  here  appeared  as  a  publication, 
the  latter  was  substituted  for  the  original  text.  A  rough  draft  of  this 
report  was  handed  out  to  Seiler  lab  personnel  on  2  Nov.  1984. 


II.  A.  Low-frequency  noise  in  ring  laser  gyros 
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Abstract 


•  i  ,*  v.-  of  rin.ilaser  beat  frequencies  Af(t)  consist  mainly 
no : se  }'  l/-  noi3e.  Consequently,  experimental  results 
■s  -  f  trie  Fourier  frequency  f, 


•  -a  n ■  ,e ,  h_,  is  usually  somewhat  above  the  Jagnepai n- 
.•  try  wi  iely  between  devices.  The  most  important  1/f 
c  t  - ring  laser  performance.  All  analysed  noise 
:  f  tie  magnitudes  of  the  noise  involved. 

•  ?  :  :  rciue.ooy  noise  in  oscillators 


>  its  ragged  appearance  on  TV  screens,  scopes, 
spectral  densities  are  plotted  versus  Fourier 
p  <es  (6 -functions)  are  essentially  absent  in 
r.  the  evaluation  of  spectra  versus  frequency 
-iging  time  r  ,  proper  account  is  taken  of  the 
i  tr.  arbitrary  frequency  dependence),  discrete 
individual  power  spectral  density  components 
dividual  erpors  A3/3  =  1.  Averaging  of  m 

.  i  Ad/3  =  m .  In  many  published  spectra, 
11  3impiy  be  fitted  by  smooth  spectra  if 
itely  obeyed  even  in  small  sampling  theory1’ “, 
'-sample  variances  (Allan  variances). 


:  .  1  »*  •  r j  mixes  no  exception.  tie  will  restrict  our  discussion 

i  :  -  ■  s  wi  t  n  tineir  usual  extreme  spectral  purity,  such  that  the 

:  .  :  .  •  •  _ •  e  s'lr.-es  in  the  lasing  system  rarely  leads  to  excursions 

•  .  .  a  :  "  here,  f ,  is  the  average  oscillation  frequency  for  which 

i.  •  -  ij-  :  *\xen.  In  t.ne  particular  case  of  a  ringlaser,  .two  or  more 
-i"  -<:  *ei  in  tne  same  cavity,  and  the  difference  frequency  Af  is 

•  e  -  i  n  •  •  :•»  3-urce  producing  correlated  effects  on  the  oscillations  may 

i  p  p  e  i r  .  n  Af.  T-.is  is  tne  major  reason  ring  laser  outputs  show  frequency  noise  which  is 
ers  -f  ;  if.  tile  c  e  1  -  w  even  tne  best  maser  clocks.  Indeed  the  irreducible  quantum  noise 
:  tut  •  r.  •  its  seen  reicoed  in  good  ring  Lasers  (gyros)  over  many  hours  of  observation 

c n  s  lute  i  remar <t3.e  feat. 

In  <eep:  r.g  wi  tn  t.ne  remarks  above  and  with  common  practice  in  oscilLator  work  we  propose* 


;  r  v  e  i 
l  ien  cy 


il  power 


1  aw  ser i es 


for  the  (one-sided)  power  spectral  density  3. 


or  the 


frequency  Af(t)  between  two  (or  more)  optical  oscillations  versus' 


r  h_2f-2 


i-  h2f 


:n  an  equivalent  Allan  variance3 

t2(t)  =  (h0/2)T~'  i-  (2  ln2  h_,  )t°  4-  (4ir2n_2/b ) T  +  1  4-  .  . 


A  given  power  law  in  3  corresponds  to  a  power  law  in 


Vote  in  particular  that  a  noise 
iuces  a  minimum  Allan  variance 


cotral  density  proportional  to  h_1  (flicker  noise)  produces  a  minimum  Allan  variance 
•3U3  r  which  i3  called  flicker  floor,  a  main  characteristic  in  oscillators  which  is  always 
served,  wi  tn  the  possible  exception  of  hydrogen  masers-'. 

To  avoid  possible  confusion  it  may  be  noted  here  that  the  noise  i'.scussei  ’  n  thi3  paper 
strictly  frequency  noise  which  i3  measured  by  counting  the  output  over  a  given  time 
terval,  or  alternately  by  comparing  the  frequency  directly  wi  tn  a  frequency  reference. 
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m  a-  *-*>  c* 


:  n  i 


Any  imp-' 
*  nclaiei 


noise  is  disregarded  in  practise,  altnoug.n  in  principle  it  would  nave  to  be 
t.norougn  discussion  of  oscillator  noi3e. 


In  tr.e  following  we  present  data  obtained  from  optical  oscillators  over  t.ne  last  two 
decades.  Tney  are  shown  either  as  a  frequency  deviation  versus  tine  from  a  reference,  as  a 
power  spectral  density  of  the  beat  frequency,  or  as  an  Allan  variance  of  the  beat  frequency 
versus  averaging  time.  T.oey  fora  the  empirical  basis  for  the  classi  f  ’  cati  on  scheme  above. 
We  t.nen  provide  some  thoughts  on  theoretical  explanations  of  the  observed  noise.  Finally, 
projected  data  on  noise  of  large  rings  are  included. 


Q 


Experimental  evidence  of  laser  frequency  fluctuations 


Some  preliminary  i 
di cussed  in  detail  in 
fluctuations  by  Jaseja 


nformation 
ref.  S. 
et  al.  3  i  s 


is  given  in  ref.  7.  Time-domain  representations  are 
One  of  the  earliest  documentations  of  optical  frequency 
gi  ven  in  f i gure  1 . 
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Figure  1 .  Fluctuation  of  frequency  difference 
between  two  lasers  at  \=  1.1 5nO 
observed  by  Jaseja  et  al.  in  1963- 


Figure  2.  Fluctuation  of  frequency 
difference  between  two 
counterci rculati ng  modes  at 
\  =  633  nn  in  a  ring  laser 
gyro  (1932). 


The  beat  note  stayed  audible  over  several  minutes 
in  contrast  to  this  is  figure  2‘  obtained  on  a 
observed  not  to  exceed  +  0.7  Hz  over  2  1/2  days, 
a  sampling  period  of  100  s  to  avoid  significant 
deviation  was  o  =  0.15  Hz  over  the  whole  period. 


,  i.e.  it  was  of  the  order  of  kilohertz, 
ring  laser  in  1932,  where  the  beat  was 
The  beat  frequency  has  to  be  oounted  ever 
quantization  errors.  The  resulting  rms- 


Turniag  to  power  spectral  measurements:  3  versus  power  input  P  into  a  mode  was  studied 
by  Kuvatova1 3  in  1976  where  a  P- 1  dependence  was  clearly  observed  over  more  than  a  decade  of 
power,  at  Fourier  frequencies  where  the  laser  noise  was  white,  fig.  3- 

As  early  as  1967,  Siegman  had  already  succeeded  in  separating  the  noise  in  a  white  part 

and  a  f-  -part,  with  a  transition  frequency  of  about  1  kHz,  using  two  separate  lasers  under 

very  careful  conditions,  fig.  4. 

The  following  figures  (5  through  10)  were  taken  from  data  on  a  commercial  ringlaser 
between  I960  and  1932. 

Figure  5  shows  a  spectrum  with  a  clear  1  ft  part  over  one  decade.  The  transition 
frequency  is  at  5  millihertz.  The  open  circles  show  the  effect  of  quantization  noise  as  an 
additional  noise  power  spectral  density  for  this  run  with  sanpling  period  10  s.  Correcting 
f;r  it  gave  the  white  ncise  as  indicated.  A  different  run  with  ICO  s  sampling  interval  shew s 
e  -illusively  Iff  r.cise,  fig.  6.  The  error  tars  are  the  statistical  (ras  1  errors  cf  clusters 

ICO  points  '±  1?S).  In  fig.  7 ,  a  long  run  Cover  9's  cays’1  shows  f”  ,  and  white  noise 

exclusively.  The  noise  nas  a  f“  behavior  over  an  impressive  five  decades  of  power  density, 
with  no  discernible  deviation  from  equation  (1)  except  possibly  below  3  mHz.  Figure  3  gives 
a  quantitative  account  of  the  relative  spectral  purity  of  ringlaser  outputs  (S^f/f  ) 
compared  to  data  or.  hydrogen  masers  given  by  Vessot  et  al.'^  in  1  977.  Note  that  the  stalls 
f  the  art  in  both  hydrogen  masers  and  ringla3ers  has  considerably  improved  since;  the 
iicker  floors  are  estimated  to  be  about  two  orders  of  magnitude  lower  in  both  types  of 
evices.  One  of  tne  better  results  achieved  in  ringlasers  is  given  in  figure  9,  where  white 
■•Ise  prevails  at  frequencies  above  IOOmHz,  with  a  spectral  density  of  *  hQ  “ 
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9x10-^  Hz-/ Hz;  again,  me  overall  spectrum  can  very  well  be  represented  by  equation  .1  )■ 
iJote  ai30  that  the  magnitudes  are  up  to  a  factor  600  different  from  those  given  in  figure  7. 

A  representative  Allan, variance  versus  t  is  given  in  figure  13.  3ince  tne  transformation 
between  3  and  is  known'',  the  curve  <r^  versus  t  can  be  quantitatively  predicted  from  power 
spectral  ^“nsities,  and  vice  versa.  The  flicker  floor  from  figure  7  corresponds  to  <r,/f0  = 
4.6  x  ID-'7,  whereas  the  flicker  floor  from  figure  9  is  <rA/f0  =  5.3  x  1  3”  . 

Tneoretical  support  of  equation  '1) 


White  noise  term  h0 

Literature  abounds  at  this  time  in  calculations  on  white  noise1 It  is  calculated  from 
quantum  noise  considerations.  According  to  arguments  on  quantum  fluctuations  of  a  system 
far  from  threshold,  with  complete  inversion,  at  Fourier  frequencies  not  exceeding  f  =  f^/3 
(3  =  quality  factor  of  passive  optical  cavity),  at  quantum-dominated  oscillation  frequencies 
f.j  >>  k'T./h,  the  'one-sided)  power  spectral  density  of  one  oscillation  is  ■'*  * 


Pin  W) 


101  102  103 

f  (Hz) 


104 


Figure  3<  Power  spectral  density  of  a 

ringlaser  output  at  \=  1.15  urn. 


Figure  5-  Hun  with  sampling  interval  13  s, 
showing  quantisation  noi3e  at 
high  frequencies  [open  circles). 
Tne  line  f°  [white  noise)  i3 
obtained  after  correction  for 
tne  quantisation  nois^.  At 
lower  frequencies,  f~ '  noi3e 
is  apparent.  Hote  the  size  of 
the  errors  of  tne  data. 


Figure  4.  Power  spectral  density  of  the  beat 
frequency  between  two  lasers  at 
\  =  633  nm. 


f  (Hz) 

Figure  6.  Run  with  sampling  interval  100  s. 

Pure  1/f  noise  prevails,  well 
documented  via  the  error  bars  of 
the  points  which  are  averages  of 
130  neighboring  individual  power 
spectral  densities. 


*4  /  SPie  Vo!  487  Physics  of  Optical  Ping  Gyros  11984) 


>v-  STS'*  5  3  2  SAt(Hz2/„,) 


SA<(hT) 


VESSOT  et  al  (  1  977) 


HeNe  RING  LASER 


--soIbilger  4  sayeh' 


DORSCHNER  et  ai 
(1980) 


10‘3  10~2 


Power  spectral  density  of  a  9  ' /2  Figure  3. 
day  run.  The  bottom  line  is  the 
wnJ te  noise  after  subtraction  of 
quantisation  noise. 


r4  ,<H2)  io-3 


Comparison  of  noise  spectral 
densities  of  relative  frequency 
fluctuations  of  a  state-of-tne 
art  hydrogen  maser  in  1977  ani 
a  ring  laser  in  1930. 


gure  9-  Power  spectral  density  of  a  Figure  10.  Allan  variance  of  two  runs  on  the 

composite  5-day  run.  Equation  same  ringlaser.  The  broad  minimum 

(1)  fits  very  well  into  the  indicates  presence  of  1 /f  noise, 

data.  Quantization  noise  has 
been  removed  (open  circles). 

SAf  =  -  h0,  13 

=  Planck'3  constant  =  6.63  x  10"^  VAs^. 

difference  frequencies  large  compared  to  lock-in  frequencies,  the  quantum  noise  of  t.-.e 
o  (or  more)  oscillations  leading  to  the  beat  frequency  are  considered  uncorrelated ,  at 
a3t  in  second  order,  so  that  the  power  spectral  densities  simply  add.  For  example  for  a 
ur-mode  ring  laser  the  observed  beat  frequency  is  expected  to  have  a  power  spectral 
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density  3Af.  f0Ur  aoj-^3  =  4-  3^r-,  if  quality  factor  and  power  input  are  assumed  to  be  the 
same  for  each ‘mode.  Published-  comparisons  of  theory  and,  experiment  are  given  by  Dorschner 
et  al'  ,  'oy  Hammons  and  Ashby'-’,  and  by  3ilger  and  3ayeh  .  The  quality  factor  of  a  passive 
cavity  can  be  determined  with  reasonable  accuracy  (error  <  10<),  and  the  power  lost  in  the 
cavity  can  be  estimated  wi  fn  about  the  sate  accuracy.  In  references  7  and  14  the  measured 
wnite  noise  level  wa3  about  5 0<  larger  than  the  one  calculated  through  equation  (3).  The 
difference  is  considered  larger,  than  the  measurement  errors;  incomplete  conversion  in  the 
lasing  levels  iiy  and  has  been  suggested.  A  ratio  -‘Ip/^l  =  3  would  account  for  the 
observed  difference.  More  precise  comparisons  of  white  noise  observations  *wi th  calculations 
are  clearly  important. 

As  things  stand  at  this  time,  we  may  argue  that  good  ring  lasers  operate  at  "high 
frequencies"  within  a  factor  two  of  the  white  noise  given  by  equation  (3).  This  is,  in 
itself,  a  remarkable  statement,  as  very  few  macroscopic  phenomena  have  been  shown  to  possess 
quantum-noise  limited  behavior  (superconducting  devices  is  another  example).  In  designs  of 
future  rdnglaser  systems  we  may  indeed  use  equation  (3)  to  predict  the  "high  frequency" 
lower  limit,  which  actually  may  be  correct  down  to  frequencies  of  the  order  of  a  reciprocal 
day,  see  ref.  3.  We  do  not  anticipate  any  dramatic  reduction  of  this  level. 

Flicker  noise  term  h_^ 

As  opposed  to  the  term  h0,  the  magnitude  of  h_1  ha3  not  been  assessed  adequately  as  to 
the  dependence  on  parameters  of  the  laser  system.  Low  frequency  fluctuations  of  this  kind 
are  sometimes  called  "bias  instability,"  from  the  misguided  notion  that  flicker  noise  may  be 
separable  from  bias  fluctuations. 

A  quick  comparison  with  ma3er  clocxs  (figure  3)  would  have  us  predict  flicker  floors 
vastly  in  excess  of  the  ones  actually  observed  in  ringlasers  (see  refs.  4  and  7);  no  doubt 
some  cancelling  of  correlated  flicker  noi3e  is  responsible. 

In  an  attempt  to  look  for  reasons  for  flicker  noise,  quartz  oscillators  may  serve  as  an 
example,  lag.nepain  et  al  °  investigated  a  large  variety  of  quartz  oscillators,  at  resonance 
frequencies  1,  MHz  to  25  MHz,  and  at  temperatures  from  1  X  to  300  K.  They  found  the 
remarkable  i--*  dependence  given  in  figure  11 


to5  10®  10 7 

Q 


Figure  11.  1 / f  component  of  power  spectral 

density  observed  at  f  =  1Hz  in 
quartz  oscillators  versus  the 
qualify  factor  of  the  quartz 
crystals. 


r(») 


Figure  12.  estimated  Allan  standard  devia¬ 
tion  of  a  large  ring  (58  nr), 
and  anticipated  sizes  of  varia¬ 
tions  in  the  earth  rate  at  time 
scales  of  10^3  (1/4h)  to  10^3 
(30  y).  T  =  tidal  effects,  1  = 
various  lunar  effects,  3,  C  = 
Suler  reap.  Chandler  wobbles. 


wr.cn  extends  over  six  decades  of  3,  suggesting  a  fundamental  role  of  the  qualify  factor  in 
the  'if  noise.  The  theoretical  approach  to  explain  this  dependence  rests  on  the 
'.ypo  thesis'  '  that  phonon  relaxation  times  t  in  a  lattice  have  a  1/f  power  spectral  density 
by  themselves,  i.e.  3^7-  f“  .  3uf  the  quality  factor  in  a  quartz  cavity  is  related  to 
*he  losses  via  viz,  y 
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Taking  the  power  spectral  densities  on  both  sides  produces  indeed  a  Q  '-law.  jagnepa:  r.  et 
al.  moreover  propose  that  the  proportionality  factor  be  in  t.ne  neighbornood  of  one,  i.e. 


3Af/f  =  Al~4/f 


with  A  =  1 . 


Wnetner  tr.is  hypothesis  can  be  successfully  transplanted  into  laser  systems  reaains  t o  be 
seen.  An  investigation  of  the  following  aspects  of  observed  1  .'f  noise  is  desirable: 

1.  Dependence  on  the  losses,  or  quality  factor  [is  Q  '  apparent?  i 

2.  Dependence  on  beam  power  (whicrs  way  to  enter  the  losses0) 

3-  Dependence  on  bias  [is  1/1“  largely  compensated?) 

4.  Dependence  on  cavity  construction  (temperature  coefficient,  elast'c  properties,  mass?) 


At  this  time,  comparing  observed  1  /  f  noises  in  figures  7  and  3  tentatively  wi  tn  the 
Dagnepai  n-'Jebersfeld  law 

°Af  =  fo  *  f  ’  b 

we  note  that  equation  (6)  predicts  lower  1 /f  noise  by  about  an  order  of  magnitude. 

Indirect  evidence  of  the  crucial  role  of  tne  laser  cavity  quality  factor  is  the  fact  t.nat 
commercial  ringlasers  with  lossy  elements  in  the  cavity  [Faraday  rotators,  e.g.'  a-- 
i  nferi or  in  1/f  noise  ("bias  stability":  to  ringlasers  without  intracavity  elements. 

Random  walk  term  h_ 2 

T.ais  term  is  known  to  be  quite  apparatus-dependent'’,  as  opposed  to  one  idea  that  the 
atom: c  system  itself  may  produce  it.  As  an  example,  temperature  fluctuations  may  possess  a 
relaxation  spectrum  of  tne  type 

"at  ; 0 ' 

For  large  thermal  relaxation  times  this  spectrum  essentially  provides  a  f”“  source  wnich  may 
be  transduced  t.irough  the  temperature  sensitivity  of  the  system  into  beat  frequency  noise 
with  an  f““  -  dependence.  Quantitative  estimates  of  this  effect  depend  crucially  on  details 
of  t.ne  construction  of  the  laser  cavity. 

Summary 

Three  major  noise  phenomena  are  well  documented  in  ring  lasers:  White  noise,  flicker 
.noise,  and  random  walk  noise.  The  white  noise  is  reasonably  well  linked  to  quantum  noise  of 
the  oscillator.  The  second  type  is  tentatively  linked  to  the  quality  factor  of  the  cavity 
via  i  f  spectra  of  the  losses;  the  latter  play  a  major  role  in  most  theories  trying  to 
explain  ' -'f  noise  in  physical  systems.  In  ringlasers,  a  large  power  law  dependence  on  I 
would  lend  support  to  Handel’s1'  theory  on  1/f  noise.  The  random  walk  3  an  be  linked  1 3 
fluctuations  in  the  support,  e.g.  due  to  transduced  temperature  fluctuations. 

Conclusions  and  suggestions 

From  t.ne  foregoing  it  is  clear  that  very  high  Q  cavities  should  be  attempted,  not  only  to 
decrease  the  quantum  noi3e  but  even  more  so  to  decrease  the  1  f  .noise.  Tne  application  of 
these  ideas  to  large  ringlasers,  e.g.  to  the  53  •sr  ring  presently  in  the  design  stage  at  the 
3eiler  lab  ( AF  Academy.  CO)  may  prove  useful.  Estimated  Allan  standard  deviations  for  such 
a  ring  are  shown  in  figure  12:  At  Fourier  frequencies  of  1  mHc  down  to,  say,  0.1  uHc,  the 
dominating  role  is  not  played  by  quantum  noise  but  by  1/f  noise.  Quantum  noise  is  expected 
to  be  negligiole  in  this  range.  The  various  points  in  f f gure  12  indicate  anticipated 
magnitudes  of  variations  in  the  rotation  vector  of  the  earth  at  the  periods  indicated.  They 
may  then  become  visible  in  the  ring  output,  should  the  estimated  noise  materialise. 
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II.  A.  2  Flicker  noise  in  frequency  fluctuations 
of  lasers  (ref.  8) 

A  survey  of  some  experimental  data  on  fluctuations  at  low  Fourier 
frequencies  suggests  a  Q'4  dependence  of  the  1/f-noise  on  the  quality 
factor  Q  of  the  passive  cavity.  We  consider  this  a  major  information 
for  designers  of  ringlaser  gyros  about  long-term  noise  of  such  a  device. 
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Measured  power  density  spectra  S^f  of  frequency  fluctu¬ 
ations  6f(t)  in  ringlasers  show  S5f  «  1/v  at  low  Fourier 
frequencies  v,  while  quantum  noise  prevails  at  higher  fre- 

2  4 

quencies.  S(5fr=(Af0/Q  )  ( 1  / v)  has  been  found  for  the  depend¬ 
ence  of  this  1/v  noise  (Flicker  noise)  on  the  quality  factor 
0,  with  A “4  (fQ=laser  frequency).  The  Q-4  dependence  is 
readily  explained  by  loss  (or  gain)  fluctuations,  using  a 
Van  der  Pol  oscillator. 


PACS  numbers:  05.40.+j,  06. 30. Ft,  42. 60. Da,  78.90.+t 
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Noise  in  lasers  has  been  a  subject  of  research  for  the  last  two  decades 


[1].  The  simple  equation  is 
3 


'5f(v)  - 


hf 

o 

— T~ 

Q  P 


(1) 


(h  is  Planck's  constant,  f0=u0/2ir  is  the  resonance  frequency,  0  is  the 
quality  factor  of  the  passive  cavity,  and  P  is  the  power  loss  per  mode) 
for  the  one-sided  power  spectral  density  (PSD)  S6fr  per  mode  versus  the 
Fourier  frequency  v.  Here  the  frequency  fluctuation  6f(t)  due  to  quantum 
noise  gives  an  approximate  white-noise  level  which  has  been  verified  by 
experiment  [2,9].  Eq.  (1)  was  generalized  by  Haken  [3], 

Let  us  consider  a  noise-driven  Van  der  Pol  oscillator  as  a  model 
for  a  laser  oscillator  [4], 

x  +  [r  -  (g  -  yx  )]x  +  u>0  x  =  N(t)  ,  (2) 


where  x  is  the  mode  amplitude,  r  is  energy  decay  rate,  g  is  the  unsatu¬ 
rated  gain,  y  is  the  saturation  parameter,  and  N(t)  is  the  noise  source 


due  to  the  spontaneous  emission  processes. 

From  this  model  one  can  find  the  PSD  of  the  frequency  fluctuations 


6f  8tt  (g-r)  n2-n i ( g2/g i )  g(f0)  thj  • 

where  and  N2  are  the  populations  of  the  lower  and  upper  laser  levels 
respectively,  and  g£  are  the  level  degeneracies,  g(f)  is  the  tran¬ 
sition  lineshape,  and  n^  is  the  number  of  thermally  emitted  photons. 
Under  limiting  conditions  (N2  >>  Nj ,  v  <<  fQ/0,  hfQ/kT  >>  1),  Eq.  (3) 
converges  to  Eq.  (1 ) . 

To  realize  that  noises  with  Fourier  frequency  dependent  PSD,  e.g. 
Stf(v)  «  1/v  (Flicker  noise)  exist,  we  consider  a  noise  mechanism  from  a 
different  source.  Let  us  assume  that  there  exist  fluctuations  in  the 
loss,  r,  independent  of  the  existence  of  white  noise,  i.e.  for  the 
purpose  of  this  derivation  N(t)  in  Eq.  (2)  is  set  to  zero.  Now  we  try 
to  find  the  PSD  of  the  frequency  fluctations  due  to  loss  fluctuations. 

We  first  establish  an  approximate  solution  to  Eq.  (2)  up  to  first 
order  in  y=g-r  [5]: 

x ( t )  =  /(p  2  cos  (?u)ot)  ,  (4) 

2  2 

where  c  =  1  -  u  /16w0S  and  y  =  g-r.  The  resonance  frequency  a  - 
therefore  depends  on  the  loss  r,  viz. 


“  "  “o  ‘ 


(g-r) 


By  making  use  of  g=r  in  a  laser  in  steady-state  operation,  we  find  the 
fluctuation  5f  with  respect  to  the  fluctuation  5r  3S 
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if 

f 


(6) 


where  Q  =  m0/r. 

The  PSD  of  the  fractional  frequency  fluctuation  can  be  related  to 
the  PSD  of  the  fractional  loss  fluctuation  by  using  Eq.  (7)  as  follows 


1 _ 

4 

64Q 


S 


5r/r(v)  ‘ 


(7) 


The  PSD  of  6 f / f  is  related  to  the  PSD  of  <$r/r  by  the  Q"^  law.  This  has 
been  observed  experimentally  in  quartz  oscillators  over  six  decades  of 
PSD  [6].  The  proportionality  to  O'4  thus  is  independent  of  the  specific 
assumptions  on  the  type  of  loss  fluctuations. 

Using  Handel's  quantum  theory  of  l/v  noise  [7],  one  would  expect 
loss  fluctuations  to  originate  from  loss  processes  inside  the  cavity 
whose  elementary  cross  sections  of  interaction  with  the  electromagnetic 
field  fluctuate  with  a  1 / v  spectrum  [8].  In  this  case  we  expect 


S 


<5r/r^ 


(8) 


where  A  is  a  constant  depending  on  the  nature  of  the  interaction. 
Therefore  the  PSD  of  the  frequency  fluctuation  per  mode  can  then  be 
written  as 


S 5 f  (v) 


Af 


64Q 


(9) 


A  ring  laser  gyro  (RLG)  can  be  used  to  study  the  frequency  fluc¬ 


tuations  in  lasers  [9].  Fig.  (1)  shows  a  measured  PSD  of  a  four-mode 
RLG  versus  Fourier  frequency  (the  RLG  was  placed  in  a  thermostat  with 
100  uK  temperature  deviation).  The  white  noise  dominates  down  to 


5xl0"4  Hz,  Flicker  noise  ( 1 / v  noise)  is  from  6xl0‘4  Hz  to  4xl0-5  Hz,  and 
l/v2  noise  predominates  at  frequencies  less  than  4x10"^  Hz.  Generally 
it  has  been  observed  that  the  PSD  of  frequency  fluctuations  of 
ringlasers  obeys  the  following  relation  [10] 

_2  _  l 

S^f(v)  =  h_2  v  +  h_iv>  +  hQ  •  (10) 

We  have  collected  data  from  different  types  of  RLG,  including  fig. 
1  and  we  verified  that  1/v  noise  is  commonly  occurring  in  most  cases. 
The  white  noise  level,  hQ  in  Eq.  (11),  was  ascertained  [9]. 

Fig.  (2)  shows  the  measured  values  of  h.^,  per  mode,  versus  the 
quality  factor  Q.  The  least  square  fitted  line  in  Fig.  (2)  gives  a 
value  256  for  the  constant  A  in  Eq.  (10). 

The  experimental  evidence  therefore  suggests  for  the  1/v  noise 
2 

(11) 

The  1/f  noise  in  laser  systems  is  thus  to  be  considered  as  a 
fundamental  process  which  affects  the  output  frequency  of  a  laser.  1/f 
noise  in  this  model  is  linked  to  laser  loss  mechanisms  which  are  in¬ 
flicted  upon  the  photon  field.  In  this  sense,  it  is  an  even  more  basic 
noise  process  than  white  noise.  Of  course,  if  in  the  absence  of  loss 
the  quality  factor  goes  to  infinity,  there  are  no  fluctuation  effects 
and  correspondingly  there  would  not  be  any  1/f  noise.  (There  would  also 
be  no  white  noise  but  for  different  reasons).  We  note  also,  that  the 
quality  factor  enters  the  formula  for  white  noise,  Eq.  (1),  only 
indirectly  because  spontaneous  emission  is  at  the  origin  of  white  noise. 
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Figure  Captions 


Figure  1:  Typical  power  spectral  density  of  the  frequency  fluctuation 

<5 f ( t )  versus  Fourier  frequency  v.  There  are  three  distinct 

regimes  where  S*v0  (white  noise),  S-v"1  (Flicker  noise),  S«v“‘ 

The  error  bars  indicate  the  statistical  accuracy  of  the 
data.  This  noise  spectrum  was  obtained  on  ringlaser  #2  from  a 
5-day  run  with  frequency  measurements  taken  every  100  s. 

Figure  2:  Summary  of  1/v-ampl itudes  versus  cavity  quality  factor  Q 

obtained  in  three  ringlasers,  with  the  best-fitted  line  ‘ -2  . 

All  ringlasers  use  HeNe  and  operate  at  f0=  474  THz 
( x=633nm) .  The  error  bars  of  RLG#1  are  of  the  same  size  as  the 


ci rcles. 


II.  B.  Astigmatic  Gaussian  beam  in  plane  ring:  Stability, 
curvature  radius  of  beam,  location  and  size  ot  waist(s) 

The  fundamental  Gaussian  beam  is  the  basic  signal  in  optical 
circuits,  as  is  the  sinusoidal  voltage  in  low-frequency  circuits.  The 
following  section  reviews  application  of  ray  matrices  to  a  closed 
optical  circuit  (plane  quadrangle)  with  the  purpose  of  establishing  a 
i/mque  complex  curvature  q  along  the  ring.  Radii  of  curvature  and  spot 
sizes  are  derived,  as  well  as  locations  and  sizes  of  waists;  the  results 
are  applied  to  an  example  (published  by  the  Seiler  group)  of  mode 
matching.  In  the  latter  case,  perfect  alignment  is  presupposed. 
Various  checks  have  been  built  in  to  improve  the  reliability  of  the 
calculations.  The  latter  make  this  section  a  proper  base  for  a  computer 


program. 


Plane  Four-Mirror  Ring  Cavity: 


Waists  (Location,  Size),  Curvature  Radii, 


Mode  Matching 


Definitions:  Beam  travels  counterclockwise,  CCW.  Radius  of  curvature 
of  mirror  or  beam  is  called  positive  if  concave  in  the 
direction  of  the  beam. 

?<  >C 

7?<c 

Waist  in  branch  A8:  location  Sy^  (off  mirror  A),  size  w^.,  etc. 

Radius  of  curvature  of  beam  entering  at  A:  R^_ 

Radius  of  curvature  of  beam  leaving  A:  RA+,  etc. 

There  are  at  most  four  waists  in  the  ring,  and  either  zero  or  one  waist 
in  each  branch.  The  length  of  the  branch  following  A  is  (positive), 


etc.  The  angle  between  entering  and  leaving  beam  at  A  is  0A  etc.  Note 
that  the  ring  is  plane. 


Procedure:  Use  ray-matrices  throughout,  and  the  law  of  propagation  of 
the  complex  curvature  q: 


%ut  ”  ^  q-jn  +  B)/(C  q-jn  +  D)  . 


The  circuit  consists  of  straight  sections  and  mirrors.  For  oblique 
incidence,  a  curved  mirror  is  astigmatic,  i.e.  it  has  a  different  focal 
length  for  sheaths  of  rays  _i_n  the  plane  of  the  ring  (x)  and  normal  to 
the  plane  of  the  ring  (y). 


Straight  section  of  length  S 


1  S 


0  1 


Mirror  with  focal  length  f 


1  0 


-1/f  1 


One  leg,  including  first  a  straight  section,  and  then  a  mirror: 


0  \  /  1  S 


v-l/f  1/  \o  1/ 


-1/f  1-S/f /  . 


Applying  four  legs  to  a  ring,  starting  CCW  from  ST,  at  A+: 


Check  whether  determinant  is  1: 

A  B 
C  D 

This  is  a  nectary  condition  in  a  round  trip. 


=  !  =4J>-X 


Check  whether  ring  is  stable;  0  i  1  -  (A  +  D)^/4  <  1 


This  is  required  for  the  uniqueness  of  the  gaussian  eigenmode. 


Complex  curvature  q  at  A+: 


From  the  condition  that  q  is  uniquely  specified  at  any  point,  we  get 


qA+  =  (A  <1a+  +  8)/(c  ^A+  +  °!  *  or 


3A+ 


=  -(D  -  A)/2C  +  [B/C  +  (0  -  A)2/4  C2)]1/2 


(the  positive  root  has  to  be  chosen)  . 


The  root  is  imaginary,  if  the  ring  is  stable.  Proof: 


Because  of  AD  -  BD  =  1, 


C  =  (AD  -  1)/B,  we  get 


qA+  =  (A  -  D)/2C  +  ( | B I / 1 AD  -  1 1  )[-( 1  -  (A  +  D)2/4)]1/2. 

A1  so 

-1  -12 
qA+  =  (D  -  A)/2B  -  j[l  -  (A  +  D)2/4]l/2/  B  *  -  j(x/irwA+). 

Location  of  waist:  At  the  waist  the  curvature  radius  of  the  beam 

O 

becomes  infinity,  Rw  ®  ,  i.e.  qw  =  j ( ttw'  /\)% 

A 

By  adding  a  straight  section  SwA  to  A,  we  arrive  at  the  waist,  where 

2 

^SwA  =  ^A+  +  SwA  =  %  =  j  (irwA/X) 

or  Re(qA+  +  SwA)  =  0  ,  or  $wA  «  (D  -  A)/2C  =  -  Re(qA+) 

Check:  a)  if  SwA  >  SA,  there  is  no  waist  in  SA 
b)  if  SwA  <  0,  there  is  no  waist  in  SA. 

Size  of  waist:  Translation  from  A  to  the  location  of  the  waist  does  not 
change  the  imaginary  part  of  q,  therefore 

2 

*w A/X  =  Im(qA+)  =  [-B/C  -  (D  -  A)2/4C2]1/2. 

Radius  of  curvature  at  mirror  A  (before  and  after  reflection) 

From  above,  we  get  for  the  radius  just  after  reflection  at  A: 


Re(qA+  )  =  Ra+  =  (D  -  A)/2B,  or  Rft+  =  2B/(D  -  A)  =  [Re(qA+  )]_1 


This  can  be  positive  or  negative. 

Radius  of  curvature  just  before  reflection  at  A: 

Use  A  B  C  D  -  law: 

-1  _i  -1  2 

qA+  =  (1  qA_  +  0) / ( -f A  qA_  +  1),  or  qA.  =  RA-  '  j(x/lTWA-) 

-1  -1 
=  ^A+  +  f A 

or 

RA.  =  [fAL  +  (D  -  A)/2B]-1  =  Cf^1  +  Re(qA+1)]'1  =  [1/Ra+  +  l/fA]_1 

Check:  Given  waist  (size  and  location)  and  the  distance,  the  following 
is  true  for  one  branch: 

2 

clA+  +  swA  =  <lw  =  J’(*wAA) 
or 

-1  2  2 
[RA+  -  jU/TTW  A+)]_1  =  j(irWA/A)  -  S  wA. 

Equating  real  and  imaginary  parts,  we  get 


12.  ^ 

RA+  =  -[SwA  +  SwA(ttwa/x)2]  ,  or  generally  R(z)  =  ±z(l  +  zj Z  ) 

2  2  > 

WA+  =  WA  +  SwA/(1TWA/x)2  »  or  generally  w2(z)  =  (1  +  z2/zj^). 


Here,  R  and  w  are  the  curvature  radius  and  the  spot  size  at  a  distance  z 

o 

from  the  waist,  Zg  is  the  Rayleigh  range,  Zr  =  irw£/x,  with  w0  beiog  the 
waist  size. 

Spot  size  at  A: 

-1  -1  2 
Since  qA+  =  RA+  -  j  ( X/ttwa+)  ,  we  get 

wA+  =  [  (  X/tt)  /  ( -  Im(qA+  ))]1//2  =  wA_ 

Waists  and  radii  in  other  branches: 


Two  procedures  can  be  applied 

1)  Recalculate  A6CD  matrix  starting  CCW  off  B,  then  C,  then  D.  All 
the  equations  above  reapply,  after  all  letters  A  are  changed 
into  B,  C,  D. 

2)  Calculate  qB+  by  transforming  qA+  through  one  leg  with  matrix 

(  ')  - 

Vi/fB  i -sA/f  B 

-l 

%+  =  (qA+  +  SA)/(-fB  qA+  +  1  -  sA/fB)'  find  Re(qB+)>  Im(qB+) , 

-1  -1 

Re(qB+  ),  Im(qB+  ), 


Application:  Mode  matching  paper  by  Seiler  group 
(Appl .  Opt.  22,  pp.  2487-2491,  1933,  ref.  9) 


Given  an  almost  square,  presumably  plane,  ring  with  the  four  sides 
SA  =  78.90  cm  (use  "cm"  to  strike  a  compromise 

Sg  =  78.74  cm  in  resolution  between  S  and  w) 

$c  =  78.58  cm 

SQ  =  78.74  cm  (use  X  =  632.8  x  10'7  cm  =  HeNe  wavelength  in  air) 


Assume  9  a  =  90°  exactly. 


Find  the  rest  of  the  geometry 


v  ?;  ■.  ft 


j  1  *?  * 


7.3. 


7  '■/  C  «• 


with  5  = 


2  2 
<Sft  *  S0)l/2 


7-i  5V'  c 


=  111.468  cm,  calculate  0C  through  cos-law 


2  2 


5^  =  Sg  +  S(~  -  2  Sg  Sc  cos  9c  »  which  gives 


9C  *  90.233° 


Now,  9g 


tan-1  SD/SA  +  cos’1  [(Sr  -  Sg  -  52)/(-2Sg$)], 


=  X  *1 


which  gives 


08  =  89.767°. 

Since  it's  a  quadrangle,  6q  =  360°  -  eA  -  eB  -  9C  =  89.999°. 

The  errors  appear  to  be  aS  ~  0.01  cm,  and  the  sides  '  100  cm, 

therefore,  use  rounded  values  eA  =  90.00° 

eB  =  89.77° 

0C  =  90.23° 

0q  =  90.00°  (Trapeze) 

Mirror  radii  (nominal ly) : 

rA=rc  =  600  cm  fAx=(rA/2)  cos( 9A/2 ) =2 1 2 . 1 32  cm,  fCx= ( rc/2 )cos( 9c/2 ) =21 1 . 706  cm 
r3  =  rD  =  fAy=(r'A^2)/cos(eA''2)=424-264  cm,  fCy  =  (  rc/2)/cos(0c/2)=425.118  cm 

Round-trip  matrix  at  A: 


i-(sc+sD)/fc 


(SA+SB )  ( 1-(Sc+Sq ) /fc )+sc+sd 


'A  B 


-(l-(Sc+SD)/fc)/fA-l/fc  (l-(Sc+SD)/fc)( i -(sA+sB )/fA)- (SA+SB ) If  aJ  V  d 


In  x-plane 


in  y-plane 


V.x 


+0.2569  +197.817 


■0.00593  -0.67713 


’A+,y 


'+0.6299  +256.623 


v-0. 003837  +0.02432 


Checks:  |MA+>xj  =  1=1  MA+jy,' 


Sx  =  1  -  ((A  +  D)2/4)  =  0.96 


Sy  =  0.89. 


Stable,  however,  close  to  limit  given  by  4  flat  mirrors:  (S  =  1), 


Calculate  qA+^x,  °lA+,y  *  anc*  their  reciprocals 


qA+jX  =  -78.694  +  j  164.743 
qA;Jx  =  -0.002361  -  0.004942 


qA+,y  =  -78.917  +  j  246.277 
q A;Jy=  -0.001180  -  j  0.003682, 


from  which  the  locations  and  the  sites  of  the  waists  are: 


S wA,x  =  -Re(qA+,x)  =  78*694  cm  SwA,y  =  +78.917  cm 

wA,x  =  C ( X/u) Im( qA+  ^ x ) 3 1 =  0.05761  cm  w A>y  =  0.07043  cm. 

Curvature  at  A: 

-1  -1  . 

rA_>x  =  +424.95  cm  =  [fAjX  +  Re(qA+jX) ]_i ,  RA-,y  =  +  849.59  cm 

-1  , 

RA+,x  =  ~423 .58  cm  =  [Re  (qA+x)]"1  ,  RA+,y  =  -847.47  cm. 


Spot  at  A 


WA+  x  =  0*^6384  cm  ,  w^+  y  =  0.07396  cm. 

Round-trip  at  C 

(since  the  mirrors  at  8  and  D  are  flat,  there  can  be  two  waists  only) 


Exchange  Sc  *  SA  ,  SB  *  Sg,  fA  *  fc,  and  feed  in  numbers 


+0.2569  +198.052 


'+0.6284  +256.506 


MC+,x  "  [  j  MC+,y=1 

-0.005927  -0.67712/  \-0. 003835  +0.02582/ 


Checks:  |MC+>X|  =  1 


!MC+,yl  =  1 


Sx  =  0.96  ,  Sy  =  0.89  as  above. 


qc+jX  =  -78.786  +  j  164.941 
qC+1x=  -0*002358  -  J  0.004936 


qC+,y  =  -78*563  +  J  246.390 
qC+!y=  -0*001175  -  J  0.003684 


3wC,x  =  78.79  cm 

w =  0.05764  cm 

Rc-X  =  +  422.73  cm 
RC+,x  =  "  424.10  cm 

WC+,X  =  0.06388  cm 

Note  that  >  Sq,  i.e.  there  i 

branch  0:  SwC>x  +  SwD>x  =  SwC>x  - 

On  the  other  hand,  $wC,y  <  $C» 


SwCjy  =  +78.56  cm 

W£#y  =  0.07045  cm 

R c — y  =+849.18  cm 

Rg+ #y  =  -851 .30  cm 
WQ+>y  =  0.07394  cm 

s  no  waist  in  the  branch  C,  but  in 
Sc  =  (78.79  -  78.58)  cm  =  0.21  cm. 
it  actually  right  at  the  mirror  0. 


Generally,  the  tangential  and  the  sagittal  waists  are  not  at  the  same 
1 ocations. 


Summary  for  x-plane 
(all  data  in  cm)  ^  + 


Waist  in  branch  A:  Sw^  =  78.92  cm  (about  at  mirror  B),  wA  =  0.0704  cm 
Waist  in  branch  C:  SwC  =  78.56  cm  (about  at  mirror  D),  wc  =  0.0704  cm 
at  mirror  A:  R^*=  +850  cm,  R^+  =  -847  cm,  wa+=wA-  =  0.0740  cm 

Rr_  =  +849  cm,  Rr+  =-851  cm,  WC+=WC-  *  0.0739  cm 


at  mirror  C: 


ZyOL  =  4-34  cm  wOL  =  0-03014  cm  x  =  126.76  cm  ZWqr-Z^=0.29  cm 
=  ^WOL  +  x  "  (^WOR  ~  2^)  =  ^-34  cm  +  126.76  cm  -  (0.29  cm)  =  130.81  cm 
^■WORx  =  x  +  Z^ql  +  126.76  cm  +  4.34  cm  =  131.10  cm 
WqR  =  0.05764  cm 

(ZML  1S  evaluated  [p.  248 8],  ZM  is  measured,  ZW0R  -  Zf1  is  evaluated) 

To  find  Z|_  =  Lens  spot,  where  the  spot  sizes,  extrapolated  from  either 
side,  are  equal  (note  that  Zw0L,  Zm,  ZWOr  are  fixed) 

2  2 

WL  =  W0L+  (x/7tW0l)2(Zl  -  ZWQL)2  =  W0L2  +  (X/ttWol)2z[Z  ,  (Zl  =  ZL  -  ZWQL) 
and 

WL2  =  W0R2  +  (x/,frW0R)2(ZW0R  ‘  ZL)2  =  W0R2  +  (xZitW0r)2  (x  ~  ZL)2, 


from  which,  by  subtraction 


W0L2  -  W0R2  =  (X/tt)2c(x  -  ZjJ2/W0R2  -  (Z‘2/W0L2)], 
or 

W0L2  W0R2(W0L2  "  W0R2)  =  U/*)2[W0L2(x  -  Z|J2  -  WQR2  Z^2] 
or 

<wol2-wor2)  zl2“2xWol2  zl+^2  w0l2+U/*)2(w0l  W0R)2  (Wor2  -w0l)]=o, 

which  is  a  quadratic  equation  for  Z^. 

(Note:  ZL  =  "Z"  in  manuscript;  two  typos  in  eq.  2) 

Solution  of  quadratic  equation  =  position  of  lens  with  respect  to  laser 
waist. 

ZL  =  ZL  ’  zW0L  =  '  x  w0L2/(W0R2  ‘  w0l2)  +  ^x2  W0L2/(W0R2  “  W0I_2) 

+  ("/x)2  W0L2  W0r2  +  x2  W0l4/(W0r2  -  W0l2)2]1Z2 

Horizontal  case: 

With  x  =  126.76  cm,  (note:  all  in  cm) 

WoRx  =  0.05764  cm, 

W0L  =  0.03014  cm, 


,»  v  V 
%  *,  % 


■Jji  mV.  m  ;.j 


\  =  632.8  x  10' '  cm. 


We  9et  Zl>x  =  ZLx  =  ZW0L  =  77.82  cm 
or  ZLx  =  77.82  +  4.34  =  82.16  cm 

Vertical  case: 

From  above,  the  vertical  waist  is  located  =  78.92  cm  CCW  from  A 
(mirror  rg),  which  makes  it  appear  in  branch  B  (S4),  at  S4  - 
SWA  =  (78.90  -  78.92)  cm  =  (-)  0.02  cm,  off  mirror  B  ( rj ) .  Projecting 
this  waist  to  the  left  of  mirror  B  by  0.02  cm,  we  get 

x  =  (ZM  -  0.02  -  4.34)cm  =  126.45  cm,  ZWQRjy  =  x  +  Zwql  =  130.79  cm, 

Wqry  =  0.07043  cm 

Zj^  =  96.13  cm,  or  ZLY  =  (96.13  +  4.34)  =  100.47  cm. 

I 

(If  the  paper's  number  =  0-21  cm  is  used,  one  gets  Z|_'y  = 

96.11  cm  orZL>Y  =  100.45  cm,  essentially  no  difference.) 

The  cylindrical  lenses  for  the  mode  matching  are  spaced  (100.45  - 
82.16)cm  =  18.29  cm  =  Z|_jY  *  ^L,X’ 

Curvature  radii  of  matching  lenses: 

-1  ,  , 
fx  =  R‘i(from  laser)  -  R"1(from  ring,  x),  with 


45- 

R ( from  laser)  =  Z^x  [1  +  (ttW^xZ^ )Xfj=  103.96  cm. 

R( from  ring,  x)  =  -  (ZW0R>X  -  ZL>X)(1  +  ( *W0R ,x/x(Z0R ,x  "  zL,x))2)  =  "604*84  cm’ 
or  fx  =  +38.71  cm. 

fy1  =  R-^from  laser)  -  R'^(from  ring,  y),  with 

'•  r  2  j  a_ 

R(  f  rom  laser)  =  ZL  y[l  +•(  itWqLA  Z^y)  J  =  117.29  cm. 

R(from  ring,  y)  =  -( ZWQR  y  -  ZLjy)(l  +  UWQRjy  AUwOR,y  -  ZL>y)2)  =-2030.50 

cm, 

or  fy  =  +110.88  cm 


The  data  agree  with  those  given  in  Appl .  Opt.  (paper  by  Seiler  group) 


Given  the  distance  between  two  waists,  and  their  sizes,  calculate  d^, 
d2,  f  focus  length  and  location  of  matching  lens.  Since  the  curvature 
radius  of  the  beam  at  a  waist  is 


qi  •  J 


q2  ■  J 


The  matrix  between  q^  and  q2  is 

A  8  W  1  d2\/  1  0\/l  dx  \/l-d2/f  d1+d2(l-d1/f)  \ 

C  D/  \ 0  l/\-l/f  1  /  \0  1  /  \-l/f  1-di/f  ) 


d2  d  j 

_  Aqt  +  B  (1  -  — )q  1  +  dx  +  d2(l  -  — ) 

Q2  =  CqT  7  D  =  7  " 


0  1 


qi  +  ! 


1  d  l  d2  d  j 

-  j  (q  1  q2)  +  (1  -  — )q2  =  (1  -  — )q i  +  dj  +  d2(l  -  j~) 

Using  the  fact  the  qj,  q2  are  imaginary,  we  compare  real  and  imaginary 
parts 

l  ir  2  2  d  j 

Real  parts  -  y  ( - j  wx  w2)  =  d1  +  d2(l  -  — )  to  determine  f  and  d:; 


d  j  ttw2 


2.  *  .2 


Imaginary  parts  (1  -  7-)  — -  =  (1  -  ~r'1  ,  W, 


From  the  latter  equation 


*<• 


2  2 
d  iW 2  +  d  2w  j ) 


2  2 

Wi  -  w2 


or  f  - 


2 

d2Wi 


w 


1 


2 

d  lw2 

2 - 

w. 


Substituting  into  real  part 


it  2  2 

~~7  w  lw2 

X 


2  2 
VI 1  -  VI 2 

- 2 - 


d,  +  d- 


d  ,w 


1  "  d  iw  2 


2 

w  1  - 


2 

1  "  w2 

- 2 - 7 

d  2w  j  -  d  jw2 


d  1  d- 


w  4,  4  4  4  4 

or  — 7  w1w2(w1  -  w2)  =  L  ( d  2w !  -  d^)  -  d1d2(w1  -  w2)  , 

X 

(L  =  dj  +  d£) 

2  2  22^2222 

or  L[(L  -  d^v*!  -  d1w2]  -  d^L  -  d^^  -  w2)  =  ^  w1w2(w1  -  w2) 

x 

2  2  2  _  2  2  2  2  7f2  2  2  2  2(  2  2 
diLwj  -  w2J  +  dJ-l-W!  -  Lw2  -  lwL  +  Lw2]  =  —5-  w jw2 ( w j  -  w2)  -L  w:, 

X 

2  22  2  22  2  2,222 
or  (Wj  -  w2)d  L  -  2Lw  j  +  [L  Wj  +  (w2  -  Wj)(-)  w^]  =  0  , 

Which  is  identical  to  the  equation  above,  except  for  notations. 

Note:  The  material  in  the  section  above  has  not  been  submitted  for 

publication,  as  no  originality  is  claimed.  However,  we  consider  the 
inclusion  of  the  many  checks  above  a  considerable  improvement  over  the 
published  papers  in  ref.  9,  and  rather  simple  for  programming. 


II.  C.  Ray  matrix  approach  to  square,  plane  ring  with  equal  mirror 
radii:  Stability,  waist  sizes,  spot  sizes  on  mirrors,  design  for 

minimum  spot  sizes  on  mirrors  (minimum  diffraction,  smallest  mirror 
size),  design  for  circular  spot  size  on  mirrors  (near-confocal  cavity 
for  best  mode-matching). 

The  results  of  the  previous  section  are  applied  to  a  highly 
symmetric  configuration.  A  square  ring  has  certain  advantages:  It  is 
easier  to  survey  than  a  general  ring,  the  waists  are  in  the  center  of 
the  branches  (if  four  equal  mirror  radii  are  employed);  there  are  no 
"holes"  in  the  stability  vs  mirror  radius;  manufacturing  of  four  equal 
mirrors  is  easier;  due  to  the  symmetry,  the  sagittal  and  tangential 
waists  are  coincident  in  the  center  of  each  branch;  there  is  also  the 
poxfibil ity  to  produce  a  circular  spot  size  on  each  mirror:  Since  the 
injection  of  an  external  beam  will  take  place  at  a  mirror,  mode-matching 
is  particularly  facilitated  in  this  case  by  making  the  backside  of  the 
mirror  into  a  cylinder  lens.  The  beam  will  then  be  strictly  stigmatic 


outside  the  cavity. 


w  - 

"ST 

* 

d 

f 

\ 

_ c 

Ray  matrix  approach  to  square,  plane  ring  with  equal  mirror  radii: 
Stability,  waist  size,  spot  sizes  or  mirror,  design  for  minimum 
Spot  size  on  mirrors  (minimum  diffraction,  smallest  mi rror  size), 
design  for  circular  spot  size  on  mirrors  (near-confocal 
cavity  forest  mode-matching). 

f/  d  *  -Nf  Relation  between  f  and  R: 

/f 

In  plane  of  ring  (p-plane): 

^  d  fp  =  (R/2)  cos  45°  =  (R/4 )  \T 

fx  1/f  Perpendicul  ar  to  plane  of  ring 

\  ^  "/  (s-plane): 

fs  =  R/(2  cos  45°)  =  (R/2 )  f? 

One  loop,  starting  CCW  at  st: 

A  1  0\  1  d\.'l  Oy  1  d  V  1  c\,'l  d  V'  1  0  |  1 

c  0/  i -f-1!  /  I  o  l/i-rh'lo  i  i -f-1!  '  (  o  l/l-rhMo 


c  d/  \-rii/\o  i/^r-4  \o  l  f\~rh> 

with  x  =  d/f,  we  get 

A  b\  1  -  6x  +  5x2  -  x3  d(4  -  lOx  +  6x2  -  x3) 

0  D J  ^d_1(-4  x  +10x2  -  6x3  +  x4)  1  -  lOx  +  15x2  -  7x3  +  x4 

Since  the  beam  "exits"  where  it  enters, 

A  B 

=  1  =  (l-6x+5x2-x3) (l-10x+15x2-7x3+x4)-d(4-10x+6x2-x3) 

C  D  d_1(-4x+10x2-6x3+x4) 

«  x°  (1) 

+  x1  (-  10  -  6  +  16) 

+  x2  (+  15  +  60  +  5  -  40  -  40) 

+  x3  (-7  -  90  -  50  -  1  +  24  +  100  +  241 

+  x4  (  +  1  +  42  +  75  +  10  -  4  -  60  -  60  -4) 

+  x5  (.  6  -  35  -  15  +  10  +  36  +  10) 


To  insure  that  inside  the  range  0  <  x  <  4  there  is  stability  for  any  x, 


we  determine  the  extrema.  ~  ~ 

3[  (A+D )  /2  y<?x  =  0  =  (1/2)  (-  16  +  40x  -  24x^  +  4x'j)  =  2(-  4  +  lOx  -  6x^  +  xJ) 

which  gives  extrema  at  x  =  2  :  (A+D)/2  =  +  1 

x  =  2  +  \2:  ( A+D ) /2  =  -  1 

x  =  2  -  \2:  (A+D)/2  =  -  1  (marginal  stability) 

y 

i.e.  the  cavity  is  stable  everwhere  in  the  range 

0  <  x  =  d/f  <  +  4  . 

Smallest  R,  to  have  stability  in  p  and  s:  Rm-jn  =  d/  V2  (7.6m/  >2  =  5.4  m) 

Waists:  Because  of  symmetry,  they  should  be  in  the  center  between 

neighboring  mirrors. 

To  find  the  round-trip  matrix  from  center  to  center,  do  this: 


Again,  with  x  =  d/f,  we  get 

’A1  B*\  ■  1  -  8x  +  10x2  -  4x3  +  x4/2  (8  -  22x  +  17x2  -  5x3  +  x4/2)(c 

vC'  D‘.  \  d'1 ( -  4x  +  10x2  -  6x3  +  x4)  1  -  8x  +  10x2  -  4x3  +  x4/2. 

Here,  A ‘ =0 * .  Also  note,  that  A’  +  D‘=A  +  0  (Stability!) 

Check,  whether  A ' D 1  -  C ' B 1  =  1: 

(1/2)  x°  (2)  +  1/2  x5  (-8-80-80-8+2+50+102+22) 

+  1/2  x6  (+10+32+10-5-30-17) 


+  (1/2)  x1  (-16-16+32) 


+  (1/2)  x  (+20+128+20-88-80)  +  1/2  (-4-4+3+5) 

+  (1/2)  x3  (-8-160-160-8+68+220+48)  +  1/2  x8  (+  1/2  -  1/2) 

+  (1/2)  x4  (+1+64+200+64+1-20-170+132-8)  =  1  . 

Waists  at  centers  of  each  branch,  w0: 

From  previous  section: 

R  =  -28 1 / ( A ' -D ' );  Since  A'  =  D‘  at  each  center,  it  follows  that  R  -*•  ®, 
which  indicates  a  waist. 

Si  ze  of  wai st: 

f/Z  f  /  2. 

CT/A)  w  2  =  8  ■  /  [l  -  (A'+D')2/4]  =  B7  [l  -  A'2]  , 

Since  A'  +D'  =  2A 1  =A+  0,  independent  of  the  starting  point  of  the 
round  trip.  Incidentally,  the  condition  1  -  A'2  >  0  recovers  the 
stability  condition. 

At  the  center  of  each  branch  we  get  thus 

l/z 

U/X)  w2  =  (d/2) (8-22x+17x2-5x3+x4/2)/  [l  -  (l-8x+10x2-4x3+x4/2 )2J 
which  can  be  factorized,  with  the  result 

4tt  w2/  x  d=(  x-4)  ( x-2-  fz){  x-2 )  ( x-2+  Y*2)/[  ( 1/2 )  ( x-2 )  ( x-2-  '/?)(x- 2+  \*2)  )'  x(4-x)‘] 
or 

w2  =  ( xd/2 it )  \j{ 4-x)/x 

Within  the  stability  limits  0  <  x  <  4,  the  waist  Wq  decreases  steadily, 

from  w0  (x  +  0)  +  -  (flat  mirrors)  to  w 0  (x  =  4)  =  0;  See  figure  below. 

Spot  size  at  mirrors,  wM: 


(n/Xjwfj  =  8/  [l-(A+D)2/4j  =2d(4-10x+6x2-x3)/[(x-2)(x-2-  Yz)(x-Z+  )  2 )  yx( 4-x)] 
Factori zi ng 

w£j  =  2( xd/n)/  yx(4-x) 

The  spot  size  on  the  mirrors  has  a  definite  minimum  versus  x  =  d/f: 

( a/3x)[x(4-x)]_^'/2  =  0  gives  x  =  2,  such  that 

WM»  minimum  "  \J^/n  • 

The  corresponding  waist  size  at  the  center  of  each  branch  is 

w0  =  ^Ad/2  =  (1/  \1)  ijxd/ir  for  minimum  spot  size  at  mirror 

Numerical  example:  X  =  514.5  nm,  d  =  7.62  m 
— *■  wQ  =  7.90  x  10"^  m  ^  0.8  mm 
— *wM  =  1.12  x  10“2  m  —  1.1  mm; 
relation  to  full  width  at  half-power  (see  Verdeyen,  p.  60) 

i - 

WpwHM  =  In  Z  w  =  1.18  w  . 

Circular  beam  cross  section  at  any  mirror  in 
square  symmetric  ring: 

First  note  that  the  mirror  spot  size  is  symmetric  with  respect  to  x  = 

2:  Put  x1  =  x  -  2,  and  get 

w2  =  w^/  (  Ad/ir)  =  [x(  l-x/4)  ]~l/2  =  (l-xl2/4)”l/2 


i .e  w(-x‘ )  =  w( +x 1 ) . 

Now,  the  ring  is  astigmatic,  i.e.  for  a  given  mirror  radius  R,  the  focal 
length  in  the  plane  of  the  ring  is  different  from  the  focal  length 
normal  to  the  pi ane. 


From  the  above  it  is  however  possible  to  choose  a  mirror  radius 
such  that  the  spot  size  at  the  mirror  is  circular;  simply  choose 
(d/fp)— 2  =  2-(d/fs) 

or 

d/(R  cos  45°/2 )-2  =  2-d/(R/2  cos  45°)  . 

This  results  in 

R  =  3  '7  d/4  *  1.06  d  .  (d  =  7.6  m  +  R  =  8.06  m) 

To  properly  mode-match  the  injected  beam,  put  a  cylinder  surface  on 
input  mirror  and  output  mirror. 

With  the  above  choice  of  R/d,  the  spot  sizes  at  the  mirrors  are  wM 

=  wM  p  =  wM  s  =  (xd/ir)1/2  x  1  .030  ,  {only  3%  above  the  minimum 

spot  size  possible) . 
and 

d/fp  =2+2/3 
d/fs  =  2  -  2/3  . 

Note  also  that  a  problem  of  instability  in  near-confocal  cavities  (where 

n 

fj  =  (d/2)  -  e,  but  f£  =  (d/2)  +  e,  see  Siegman,  problem  7  in  chpt.  8) 

does  not  appear  here:  Given  the  nominal  mirror  radii  R,  the  actual 

focal  lengths  fp  and  fs  are  quite  removed  from  f  =  d/2,  viz. 

fp  =  (d/2) (3/4)  ,  fs  =  (d/2)  (3/2) , 

due  to  astigmatism. 

That  means  that  the  tolerances  of  the  mirror  radii  have  merely  to 
be  <  25%,  to  avoid  near-confocal  instability. 


II.  D.  Injection  errors:  Response  of  ring  to  offset  and  tilt 

of  injecting  ray. 

In  this  section,  the  excursion  of  a  pencil  ray  off  the  optical  axis 
due  to  tilting  and/or  offset  of  an  injecting  ray  is  investigated.  The 
method  is  reminiscent  of  orbit  calculations  in  synchrotrons ,  except  that 
in  this  section  no  wave  picture  is  employed.  The  solution  of  the 
standard  difference  equation  for  rays  with  such  perturbations  is 
simple.  It  gives  however  some  insight  into  the  "stability"  of  a  given 
ring  against  injection  errors. 


II.  D.  Injection  errors:  Response  of  ring  to  offset  and  tilt  of 

injecting  ray. 


(For  analytical  background  see  Verdeyen,  Chap.  2) 


» 

Assume  an  initial  offset  rQ  =  a,  and  an  initial  slope  rQ=m  at  ST. 

The  offset  r  at  any  point  (off  optical  axis)  on  the  ring  is  then  given 
through  a  difference 

r( s+2 )  -  (A+D)  r( s+1 )  +  r(s)  =  0, 

where  r(s+v)  =  offset  at  ST  after  s  +  V  round  trips.  The  offset 

anywhere  in  the  ring  is  obtained  by  choosing  a  non-integer  value  for  s. 

Solution:  r(s)  =  rmax  sin  ( s  6  *«) 

where  rmax  is  the  maximum  offset  in  the  ring, 


cos  9  =  (1/2)  (A+D)  ,  and 


tan  a  =  [1- (A+D )2/4]1/2/[ ( A-D )/2+Br^/r0] . 


From  r(s  =  0)  =  r0  =  rmax  sin  a  ,  we  determine  rmax  as 
rmax  *  (1/S)  [(I-AOJTq  *  B2r'2  +  (A-D)B  r„  r i]1/£, 

with  S  =  stability  condition  =  [1  -  (A+D)2/4]^/2  . 

(0  <  S  <  1) 

The  maximum  excursion  is  homogeneous  in  r0  and  r0y)as  it  should  be.  In 
the  limit  of  stability's  +  0,  diverges. 

The  same  difference  equation  holds  also  true  for  the  slopes,  viz. 

r' (s+2)  -  (A+D)  r' (s+1)  +  r‘ (s)  =  0 
with  the  solution 

r'  (s)  =  r'max  sin  (s5  +  <x') 
cos  0  =  (1/2)  (A+D)  ,  arid 

tan  a'  =  [l-(A+D)2/4]1/2/[(D-A)/2+Cr0/r^],  so  that 

r'max  =  (1/S)  [(l-AD)r*2  +  C2r2  +  (D-A)C  r0r*]1/2. 

Note  that  since  a  and  a1  and  the  roots  are  multivalued,  there  exist  in 
principle  solutions  for  a  and  a'  in  all  four  quadrants,  given  the 
initial  conditions  rQ  and  r^,  and  the  ring  matrix.  In  the  following  we 
will  convene  rmax  and  rrJiax  as  positive,  and  a  and  a'  shall  be  uniquely 
specified  by  choosing  them  such  that  r(0)  =  rQ,  r(l)  =  Ar0  +  Btq,  as 
well  as  r'(0)  =  Tq ,  r'(l)  =  Cr0  +  Dr,!,. 

Special  cases: 

a)  No  initial  offset,  no  initial  slope:  Proper  operation. 

rn  *  ri  =  0:  rm-y  =  r,l,,y  =  0,  the  ray  stays  on  the  optical 


b)  No  offset,  but  slope: 


r0  =  0,  r'  =  m:  rmax  =  (8/S)m,  rj,ax  =  [(-BC)1/2/S]m 

a  =  0  or  IT  ,  a1  =  tan"*  ( 25/  (D-A) ) . 

c)  Offset,  but  no  slope 

ro  a>  ro  ^ •  r max  =  t(~BC)  ^2/S]a  r(nax  =  (C/S)a 

a  =  tan"^  ( 2 S/ ( A-D ) ) ,  a1  =  0  or  u  . 

Application  to  square,  symmetric  ring  with  four  equal  mirrors 


Starting  off  one  mirror,  the  matrix  is,  with  x  =  d/f: 

A  B  \  /l-6x+5x2-x3  d(4-10x+6x2-x3)  \ 

C  ^  d“*  ( -4x+10x2-6x3+x4)  1 -10x+1 5x2-7x3+x4 

f  1-6x+5x2-x3  d(2-x)(x-2+  f7)(x-2-  /?)  \ 

(x/d)(2-x)(x-2-  l?)(x-2+  V2*)  l-10x+15x2-7x3+x4  / 


where  (B/d)x  -  -  dC, 

S  =  (1/2)  (x-2)  (x-2-  }/2)(x-2+  VT)  •  [x(4-x)D1/2  »  [1  -  (A+D)2/4]1/2 

1-AD  =  -BC  =  x(2-x)2(  x-2-  U )2  (x-2 +  vT)2 
A  -  D  =  (x/d)B  =  -dC  *  x(2-x)  (x-2-  V?)  (x-2+ y7 ) 


and,  as  before,  AD  -  BC  =  1,  with  the  stability  limits 

0  4X44. 

For  the  special  cases  above 


=  0,  Tq  =  m,  we  get  a  =  0  or 


rn  =  a,  r'=  0 


a  '=  tan-^  C(4-x)/x]^2, 

a  =  tan"*  (4-x)*,/2, 

a'  =  0  or  it  , 


rmax  =  2dm/[x(4-x)]1/2 


rmax  =  2m/(4-x)1/2 


rmax  "  2a/(4-x)1/2 


r  max  =  (2a/d)[x/(4-x)]1'/2 


Generally,  the  square  symmetric  ring  with  equal  mirror  radii  gives 


max 


*  2(xr2+d2r^2+dxr0r^)1/2/[x(4-x)]1/2,  tan  a  =  [x(4-x)]1/2(-r0)/(xr0+2dr‘ ) 


rmax  =  2[rQ2+(xr2/d2)  +  (r0r!!)/cl}]1/2/(4-x)1/2>  tan  a'  =  E(4-x)/x]1/2  rg/[ro+(r0/d)]. 


The  offset  at  a  given  mirror  is  given  through  the  relation 
r(s)  =  rmax  sin  (s  9  +  a  ),  and  the  slope  through 
r'(s)  =  rmax  sin  9  + 

where  now  s  has  a  fractional  value.  Suppose  the  matrix  is  evaluated 
starting  off  mirror  A  in  CCW  direction.  Then 
B 


rA  =  rmax  sin  (0  9  +  a) 

CD  rB  =  rmax  sin  (8/4  +a) 

rC  =  rmax  sin  (9/2  +a) 
rD  *  rmax  sin  (30/4  +  «). 

during  the  first  round  trip.  In  the  second  round  trip  increase  s  by  1, 

t 

etc.  Equivalent  relations  give  rA  ,  etc.  Three  examples  are  worked  out 
below.  All  results  are  verified  by  applying  the  matrix  of  one  leg  to 


the  column  vector,  i.e. 

(%«)  />  °\  Av) 

U ,,/  1/  \r J 


where  v  is  the  vth  corner. 


No  initial  offset  is  considered,  i.e.  with  rQ  =  0,  Tq  =  m,  we  have 
rmax  =  2dm/[x(4-x)]ly/2  ,  tan  a  *  0  (a  =  0  or  it) 

rmax  =  2m/(4-x)1/2  ,  tan  a'  =  E (4-x) /x]1 /2 


1.  Closed  path  after  one  round  trip:  e  =  2  it  or  cos  9  =  1=  (1/2)  (A+D) 
=  (1/2)  (2-16x+20x2-8x3+x4)  Solutions:  x  =  0,  2,  2,  4 


Only  x  =  2  is  inside  the  stability  limit:  confocal 


x  =  d/f  =  2  -*•  rmax  =  dm,  a  *  0  or  it 

*  rmax  =  a'  =  tan-1(±l)=  +  tt/4,  ±  tt/4  + 

r  -  dm  sin  ( s2^-)^  r 1  =  Hf  m  sin  (s2  it  +  3ir/4) 

( a1  =3ir/4^  to  get  =  r0  =  m  and  rB  =  (l-x)m  =  -m) 
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All  the  values  have  been  checked  directly. 

Example:  d  =  7.6  m,  m  =  tan  1'  =  0.00029 

rA  =  rc  =  0,  rB  =  -rD  =  2.21  x  10"3  meter. 

2.  Closed  path  after  two  round  trips:  9  =  tt  or  cos  9  =  -1  =  (A+D)/2; 
Solutions  x  =  2+  |  2,  2+  Y IT,  2-  l^,  2-  iT* 

>-  7 

x  =  d/f  =  2  -  |  2  (the  other  solution  corresponds  actually  to  9  =  3-rr) 


1  max 


=  t  2  dm,  «'  =  0 


rmax  ^  ^ 


r  =  2  dm  sin  ( sir) 


2-  )?'m,  a'  =  tan_1[±  (1+  ^*2)] 

=  ±  3  rr/8 ,  ±  3ir/8  +  tt 
r'  =  f?  |  2-  fT^sin  (s  n  +  5tt/8) 
(rA  =  m,  rB  *  (l-x)m  =  (  Y^-Om) 
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The  offsets  are  only  slightly  larger  than  in  the  confocal  arrangement: 
rmax  =  O.08)  Jm 

3.  Two  round  trips,  but  with  9  =  3* 
x  =  d/f  =  2  +  r T 
r' 


rmax  =  ^  dm 


max 


=  )  7  ’/§  +  ft 


a  =  0 

=  tanti/F  -1 ) 

=  ±  tt/8,  ±  tt/8  +  IT 


r  =  \  2*  dm  sin  ( s3  n) ,  r'  =  /7  V2+YTV,r7sin  (s3*  +  7  tt/8  ) 
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The  stronger  focussing  makes  again  slightly  larger  offsets:  r_ax 


1.41  dm 


The  figure  gives  an  impression  of  the  expected  maximum  excursion 
due  to  offset  {at  0,  m=0)  and  due  to  tilt  (rrtfO,  a=0):  Going  towards  the 
strong  focussing  limit  x  +  4,  the  excursions  diverge.  On  the  other 
side,  for  weak  focussing  (x  +  0),  only  a  tilted  ray  diverges.  It  is 
quite  obvious  that  from  the  point  of  view  of  minimizing  effects  of 
injection  errors,  intermediate  stability  figures  are  to  be  sought. 


II.  E.  Resonant  frequencies  of  Hermi te-Gaussi an  modes  in  a  ringlaser 

Although  the  working  modes  in  a  ring  laser  gyro  are  fundamental 
Gaussians,  a  knowledge  of  the  higher-order  eigenmodes  is  important  for 
design  purposes.  Imperfect  injection  into  a  ring  excites  higher-order 
modes  (section  1 1 . F . ) ;  the  resulting  pulling  effects  depend  on  amplitude 
and  frequency  separation  of  these  modes. 

Collins'  theory  (ref.  10)  is  applied  to  a  square-shaped  ring  with 
zero,  one,  two,  four  equally  curved  mirrors  (in  any  one  of  these  four 
cases,  the  other  mirrors  are  plane).  Simple  formulas  emerge  which  are 
useful  for  design. 


Resonant  frequencies  of  Hermi te-Gaussi an  modes  in  a  ringlaser 
This  section  presents  a  general  formula  for  the  resonance 


frequencies  of  a  given  resonator  with  negligible  apertures  and  symmetric 
focusing  elements  with  respect  to  two  normal  planes.  An  astigmatic, 
plane,  four  mirror  resonator  with  spherical  mirrors  falls  in  this 
category.  The  formula  is  derived  in  ref.  10.  A  configuration  is  sought 
where  by  proper  choice  of  mirror  radii  the  frequency  difference  of  those 
modes  to  the  fundamental  mode  i S  maximized,  for  minimum  interaction. 

The  resonator  mirrors  may  be  replaced  by  thin  lenses.  The 
equivalent  focal  length  of  a  thin  lens  depends  on  mirror  curvature 
radius  R  and  the  incident  angle  of  the  beam.  This  is  shown  in  Figure  i. 
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For  quadrangular  resonators,  the  resonance  frequencies  are  generally  ^ 
fqmn  =  (C/L)  [q  +  (m  +  0  •  5  )  ax  /  2  it  *  [n  +  0.5)<x  /2,r]  (1) 

The  parameters  are 

c  =  speed  of  light  in  vacuum  (the  resonator  operates  in  vacuum) 
l  =  optical  distance  coefficient  (round-trip  path  length) 
q  =  axial  mode  number,  m,  n  =  transverse  mode  numbers. 


a  .  arccos(.'-^«)l/2  ( u  =  *  or  y ) 
u 


Pull  =(l/2)(l/r1  -  |D0,  Q | / 1 0  | ) , 


=  (-1/2)  (l/r2  +!Dn.1jN|/|D|) 


(N  =  number  of  mirrors,  rj.,  r2  =  radii  of  curvature  of  the  reference 
surface. 


(-D  |Dn  J/|D|, 


PulN  "  '  'uo,nI 


f  i+f 2-d  -f2  0  0 

-f2  Vf3‘d2  -f2  °  •  • 


|D0  ol  = 


fj  and  d-j  =  shown  in  Figure  2-2. 
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flx  =  Rcos  45 °/2  =  (/2/4)R  (obliquity  angles  are  all  45°) 
fj  =  R/(2cos45°)  =  (/2/2)R 

The  curvature  radii  of  the  reference  surface  are  r^  ->•«>,  r 2  -*■  ®,  and 
N=2. 

The  angles  ax  and  ay  are  obtained  from  Eq.  2 
after  substituting  eq's.  (3),  (4),  and  (5)  as 

cos  a  = 
u 

Now  we  need  to  calculate  [Do  o|,  ! D 1  2 1 »  and  I D  0  2I.  These  determinants 


|Do»  ol  ID  1,2 
I D  0*  2I 


COS  a  =  1  -  L./<:Tu 
V 

Therefore  the  resonance  frequencies  are  given  as  follows 
fqmn  =(c/'-)[q+(m+0*5)arccos(l-L/2'/R)/2Tr+(n+0.5)arccos(l-L//2R)/2Tr] 0  (10) 
The  stability  condition  is  given  by 
-1  <  cos  a „  <  1 

or  (1/2)  <  R/L  <  ®  (the  stronger  focussing  tangential  plane 

dictates  the  lower  limit  for  R). 


Square  ring  with  two  diagonally  opposite  flat  mirrors  and  two  equally 
curved  mirrors  (figure  5). 
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The  equivalent  lens  wave  guide  for  the  square  ring  is  shown  in  Figure 
6.  Here,  again,  r^  =  r£  -*■  ®,  but  N=3  (two  curved  mirrors). 


The  angles  ax  and  ciy  are  obtained  by  Eq.  (2)  as 


cos  a  = 
u 


I  0  0 »  ol  l°2»  3l 
^0  *3! 


'/a 


where  I Dq s 0 1  *  I D2 , 3 1 »  and  I D0 , 3 1  are  calculated  as  follows 
|D0  ol  = 


2f  -2 1 
u 
-f 
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2  2 

=  2  ( f  .£)  -f 
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Thus , 


€n\A 


cos  a  -Z(l  -  -f)-l  ,i.e  cosa  =2(1-  and  cosa  *2(1-  -1 


The  resonance  frequencies  are 


f  =(c/L)[q+(m+0.5)arccos(l-L//2R)/ir+(n+0.5)arccos(l-L/2/2R)/Tr] 
qmn  v  ^ 


The  stability  condition  is  here 


(1/2/2)  <  R/L  >  °° 


Square  ring  with  four  equally  curved  mirrors  (see  figure  7), 


The  equivalent  lens  waveguide  for  the  square  ring  is  shown  in  Figure 
8.  This  example  shows  an  alternative  way  to  find  the  resonance 
frequencies.  It  is  assumed  that  the  total  phase  shift  due  to  one  round 
trip  around  the  resonator  must  be  2 -rrq  where  q  is  an  integer.  Therefore  we 
need  to  find  the  total  round  trip  phase  shift 

The  longitudinal  phase  shift  is  4>d=kd)J  (ref.  11)  due  to  the  distance 


du  =  L/2  between  a  lens  and  a;  the  transverse  phase  shift  is 


$  =  ~ ( u  +  0.5)  arctan  (d  /z0  ) 


where  u  =  m  or  n. 

Here,  the  Rayleigh-length  is  zQp  = 
The  total  phase  shift  is  therefore 


(L/8)[fu(4-L/4fu)/(L/4)] 


4>tr=4k£  -  (m  +  -jJarctanCf x(4- £/f 


1  /2  l 

,)/*]"  -(n  +  T)arctan[f^(4-£/f^)/£] 


1/2 


We  need  to  satisfy  =  2?rq;  therefore  the  resonance  frequency  is 


fqmn  =  (c/L )  [q+2(m+0.5  )arccos  (1-L/2/2R) /it+2(  n+0.5  )arccos  ( 1-L/4/2R ) /it]  (12) 
Pi scussion: 

The  three  types  of  square  rings  give  surprisingly  simple,  and 
similar,  results  for  the  resonance  frequencies.  With  a  parameter  6,  5  = 
1,2,4  for  rings  with  1,2,4  (equal)  mirror  curvatures,  respectively,  the 
resonance  frequencies  (10,  11,  12)  are  given  by  the  same  equation  for  all 
three  cases  as: 


f  =(cq/L )+(c5/2TrL)[(m+0. 5)arccos(l-/?L/SR)+(n+0.5)arccos( 1-L//2SR) ]  , 
qmn 

A  cavity  with  4  flat  mirrors  is  contained  here  is  as  a  special  case  with 
•5=4, R-*- 

The  spacing  of  the  fundamental  modes  (m=n=0)  is  always  as  usual 
Af  =  (c/round  trip  path  length) (Aq=1)  =  c/L  (free  spectral  range) 

This  is  true  for  any  6. 

The  arguments  of  the  arccosines  vary  with  R.  To  achieve  a  specific 
operation,  given  by  a  point  in  the  stability  diagram,  5R/t  has  to  be 
chosen  properly. 

The  two  stability  limits  are,  for  any  6, 
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and 


A 


k  =  2 
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(mirrors  with 
maximum  curvature) 


They  provide  us  also  with  limiting  frequency  differences. 
In  the  former  case 

★ 

f  =  (c/L)q  (all  m,n  modes  resonate  at  the 

c*’m*n  same  frequency  for  a  given  q). 


in  the  latter  case 


f*  m  n  *  (c/L )q  +  (c/4L)5[2(m  +  0.5)  +  (n  +  0.5)] 

such  that  all  resonances  are  either  0,  1/4,  1/2,  3/4  or  1  free  spectral 

range  apart  from  the  fundamental  mode. 

In  the  latter  case,  for  example,  a  (q,  0,  0)  mode  and  a  (q-1,  1,  0) 

mode  would  have  the  same  resonance  frequency  in  a  2  curved  mirror  ring  (6 

=  2). 

The  problem  is  therefore  to  find  a  ring  configuration  where  the 
unwanted  Hermite-Gaussian  resonances  are  as  far  away  as  possible  from  the 
desired  qOO  resonances. 

To  facilitate  the  search  for  a  ring  configuration  with  well-spaced 
resonances  we  developed  a  graph.  Fig.  9,  that  contains  the  relative 

spacing  (f  -  f  0  o)/(fn+i  o  o H2/6)  (6=number  of  equally  curved 

M  »ri  q ,  ,  q^i  ,  # 

mirrors)  versus  a  parameter  C=2/2L/sR  within  the  stability  limits  0 <C< 1 . 

It  appears  from  these  calculations,  that  not  much,  if  anything,  can 


be  gained  by  avoiding  these  special  symmetry  cases,  as  indeed  any 
frequency  distribution  pattern  can  be  identically  reproduced  by  a  square 
resonator  with  1,  2,  4  curved  mirrors;  all  what  needs  to  be  done  is  to 
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Example: 

1  curved +  3  flats  (8  =  1):  Co,,** 0.50  (R  =  5.661) 


adjust  the  parameter  (L/5R)  to  the  same  value. 

After  these  results  had  been  obtained,  a  scheme  was  developed  where 
the  (L/6R)  was  evaluated  for  optimum  discrimination  against  Hermite- 
Gaussians.  The  results  are  found,  but  not  included  in  this  report.  They 
show  that  the  problem  of  pulling  has  a  unique  solution. 


II.  F. Injection  errors:  Misalignment  and  mismatching  of  an  injected 

Gaussian  beam.  Effects  of  offset,  tilt,  mismatch  by  creating  Hermite- 


Gaussian  eigenmodes  (H.-G.'s).  Power  loss  and  pulling  of 
(fundamental)  Gaussian. 

A  real  ring  has  finite  adjustment  errors.  Since  H.-G.'s  form  a 
orthogonal  and  complete  base,  such  errors  can  be  expressed  in  an  H.-G. 
expansion.  In  real  life,  a  ring  behaves  indeed  that  way,  as  a  proper 
interpretation  of  experimental  data  in  ref.  12  shows.  The  latter  was 
done,  among  other  things,  by  us  in  ref.  13,  and  is  summarized  at  the  end 
of  this  section.  We  consider  this  section  as  a  basic  contribution  to 
design  of  a  large  ring,  and  we  intend  to  publish  it. 

The  complexity  of  our  approach  is  somewhere  between  our  ray  approach 
(section  II. 0.)  and  that  of  ref.  14,  the  latter  using  an  expansion  in 
elliptic  wave  functions.  We  treated,  separately,  three  different 
adjustment  errors:  Offset,  tilt,  mismatch.  Only  the  last  type  receives 
attention  in  the  literature,  although  it  is  actually  the  most  benign  case, 
giving  rise  to  rapidly  decaying  even-indexed  H.-G.'s.  The  other  types 
usually  have  all  H.-G.'s  in  their  expansion.  The  crucial  parameter  is  the 
ratio  of  "mi sadjustment"  to  spot  width. 

We  were  able  to  quantitatively  explain  the  data  in  ref.  12,  where 


most  H.-G.'s  were  indeed  due  to  misalignment,  contrary  to  the  claim  stated 
therein. 


II.  F.  Injection  errors:  Misalignment  and  mismatching  of  an 
injecting  Gaussian  beam.  Effects  of  offset,  tilt,  mismatch 
in  creating  Hermi te-Gaussian  eigenmodes.  Power  loss  and 
pulling  of  (fundamental)  Gaussian. 

II.  F.  1  Introduction 

In  order  to  obtain  a  detailed  description  of  the  optical  modes  in 
resonators,  the  diffration  theory  based  on  the  Huygens-Fresnel 1  s 
principle  is  used^.  The  fields  across  the  reference  surfaces  are 
extended  by  means  of  the  Fresnel-Ki rchhoff  integral  in  second  order  to 
give  the  fields  everywhere  inside  and  outside  the  resonator.  This 
integral  is  of  the  form 

ie-Jkz  -j(k/2z)[(x-x0)2+(y-y0)2] 

u(x,y ,z)  =  Yx  //  u0(x0,y0)e  dx0dy0, 

(1) 

x,y,z  =  a  rectangular  coordinate  system  where  the  z  axis  is  along  the 
direction  of  the  wave  propagation,  xQ,  yQ  =  coordinates  on  the  surface 
of  the  reference  surface  mirror,  at  zQ  =  0. 

The  parameters  are: 

u(x,y,z)  =  complex  wave  amplitude, 

u0(x0,y0)  =  complex  wave  amplitude  of  input  at  mirror, 
k  =  wave  number  =  2w/\  , 

S  =  reference  surface, 

In  order  to  find  a  solution  for  the  above  integral,  it  is  assumed 
that  the  resonator  mirrors  are  circularly  symmetric  which  basically 


leads  to  Gaussian-Laguerre  solutions,  or  they  are  rectangularly 
symmetric  which  leads  to  the  Hermite-Gaussian  solutions.  As  a  practical 
matter,  any  slight  misalignments  or  mirror  imperfections  may  inhibit 
circular  modes  (Bessel  modes), and  the  resonator  may  oscillate  in 
H.-G.'s.  Therefore,  Hermite-Gaussian  solutions  are  considered  in  this 
report  exclusively. 

The  Hermite-Gaussian  solutions  for  the  amplitude  umn  are 
C  i  j  i 

^(2V^n(4^)exP[“  ^  (*  V)  ]exp[  j  (m+n+1  )*] 

0  (2) 


The  parameters  are  defined  as: 
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w0  =  waist  size,  w  =  spot  size,  R  =  radius  of  curvature  of  beam  phase 

m  =  number  of  nodes  along  x  axis 

n  =  number  of  nodes  along  y  axis 

2 

z0  =  Rayleigh  length  =  ttwoA 
<j>  =  arctan  (z/zQ) 


*„<«>  ’  H,<Oe'{2/2 

2  n  _1 

H  (?)  =  Hermite  polynomials  of  degree  n,  H  (?)  =  (-l)T1e’  —  e"^  . 
n  n  9?n 


An  open  optical  resonator  (e.g.  Fabry-Perot  resonator)  is  capable 
of  oscillating  with  different  axial  and  transverse  modes.  As  Haus^ 


states,  the  electric  field  Emn  of  a  sperical  Hermi te-Gaussian  wave 
propagating  in  the  z-direction  is. 


3u 


En.n(x'y>z-t)  =  /C2"P:l  (xumn  ’  j 


(3) 


with  the  major  component  of  polarization  assumed  to  be  parallel  to  the 
x-axis. 

The  parameters  are: 
n  =  wave  impedance 

P  =  total  power  crossing  an  arbitrary  z-plane 
w  =  angular  frequency 
x  =  unit  vector  in  x-direction 

A 

z  =  unit  vector  in  z-direction 

Since  the  electric  field  waves  are  essentially  monochromatic  waves, 
Emn(x,y,z,t)  is  written  as 


Em„(x,y,2,t)  ■  Emn(x,y,z)e 


jut 


-V 

Emn(x,y,z)  only  will  be  used  in  the  following. 

The  z-component  of  the  electric  field  is  generally  small  with 
respect  to  the  x-component,  therefore  it  is  omitted.  The  approximated 
electric  field  is  then 

E  »  /[2nP]  u  *"jkZ  x  .  (4) 

mn  mn 


An  open  resonator  can  be  excited  by  an  external  beam.  The  external 
source  should  satisfy  the  following  conditions: 


1. 


Its  mode  should  be  transferred  into  a  fundamental  gaussian  ring 
mode  as  completely  as  possible. 


2. 


It  should  have  a  frequency  near  the  fundamental  resonance 
frequencies  of  the  resonator. 


To  fulfill  the  first  requirement,  the  resonator  has  to  be  carefully 
aligned  and  matched.  We  will  consider  misalignment  problems  in  the 
first  section  of  the  report. 

Higher-order  Gaussian  beam  modes  are  usually  undesirable  in 
resonator  applications  since  a  passive  resonant  ring  laser  gyroscope  may 
give  a  false  rotation  signal  due  to  the  presence  of  the  such  modes. 
Irrespective  of  the  amplitudes  of  such  H.-G.'s  one  tries  to  design  a 
resonator  so  that  they  appear  as  far  away  as  possible  from  the 
fundamental  mode.  But  since  there  is  always  a  finite  linewidth 
associated  with  the  resonance  frequencies,  the  tail  of  the  resonance 
does  pull  the  fundamental  mode's  frequency.  Therefore  the  maximization 
of  the  frequency  difference  between  the  fundamental  mode  (TEMqq)  and  the 
next  higher-order  mode  is  needed. 

II.  F.  2)  Mi  sal ignment 

A  resonator-source  system  can  be  misaligned  in  two  ways: 

a)  The  resonator  becomes  misaligned  internally, 

b)  The  external  source  is  not  aligned  with  the  resonator. 

The  latter  is  dicussed  throughout  the  report.  The  new  optical  axis 
of  an  internally  misaligned  resonator  can  be  found  by  applying 
variational  calculus  (minimize  pathlength).  The  external  source  can 
then  be  explained  relative  to  this  new  base  by  the  methods  outlined 
below.  It  is  assumed  that  the  resonator  fields  are  given  by  Eq.  (4). 


In  this  report,  we  consider  the  case  where  the  injected  beam  is  a 
TEMqq  mode  only.  Three  types  of  external  misalignment  are  considered. 

a)  Offset  source  beam: 

The  center  of  the  injected  beam  is  not  on  the  z-axis  (optical  axis) 
but  the  beam  is  parallel  to  the  z-axis,  with  the  same  waist  size  as 
that  of  the  resonator. 

b)  Tilted  source  beam: 

The  injected  beam  makes  a  finite  angle  with  the  z-axis  (again,  the 
waist  size  is  the  same  as  that  of  the  resonator). 

c)  Mismatched  source  beam: 

The  waist  size  of  the  injected  beam  does  not  match  the  waist  size  of 
the  resonator  beam,  but  the  beam  shall  not  be  off  the  z  axis  nor 
shall  make  an  angle  with  z-axis. 

The  above  types  of  misalignment  are  treated  in  detail  in  the 
following  subsections. 

II.  F.  2.  1)  Offset  Source  Beam 

As  it  is  shown  in  Figure  1,  the  external  laser  beam  is  a  TEMQg-beam 
parallel  to  the  z-axis  but  the  center  of  the  beam  is  off  by  e. 

The  electric  field  of  the  injected  beam  can  be  written  as 


En  =  /[2nP]  u'00e‘jkzx 


where 


U  n  n  = 


C.  f 

~2 —  exp  [- 


(x+e)  +> 


-]exp  {-  4k(x+e)2+y2]>ej*. 


It  is  easy  to  show  that  the  functions  umn(m,  n  =  0,  1,  2,  ...)  form  an 
orthogonal  and  complete  set  ^,16.  It  is  therefore  possible  to  expand 
the  injected  beam  in  terms  of  this  set,  i.e. 


E  E  A  E 

_  mn  mn 

m=o  n=o 


where 


^  ^  Emn  E0  dx  dy 

-°°  -oo  mn  v _ _ 

OO  OO  ^ 

t  r  E  E  dxdy 
/  f  mn  mn  J 


,  E  is  the  conjugate  complex  of  E 
<rm  3  mn 


To  simplify  the  analysis,  we  expand  the  injected  beam  at  the  location  of 
the  waist  of  the  proper  resonator  modes.  It  is  assumed  for  this  purpose 
that  all  H-G's  of  such  a  resonator  have  their  waists  at  the  same 
location  at  which  the  origin  of  our  coordinate  system  is  placed. 


mn  oo  oo 


/  /  (  /[2nP]  u*n  ejkz)(/[2nP]  u{,0  e'jk2)dxdyo#  m 

- — : - : — =  /  /ulu'oodxdy 

/  /  ( /[ 2nP]  u  e^kz) {/[2nP]  u  e'jkz)dxdy"“ 


at  z  =  0  we  have 


Amn  /  /  exp  [-  dxdy  (7) 


where  we  use  the  fact  that  R  -*■  »  and  $  0  as  z  -*■  0. 

The  double  integral  can  be  separated  into  two  integrals  over  x  and 
y,  respectively.  It  may  be  written  as 

(x+s)2 

2  2  2  1 


T  2  2  1 

0°  W  00  -v  /w  — 

,  V2x.  0  /2y.  y  '  0  2  ^ 

A  =  C  /  dx ( - )e  }  /  dy ( — -)e  (—5 — )  > 

mn  mn  w0  n  w0 


The  second  integral  is  unity  for  n  =  0  and  zero  otherwise1.  Thus, 

2 


mn 


(x+e) 

Cmo  /  ax  (*roOe  W°  1  for  n  =  o 


at 


for  n  i  O  . 


The  resulting  expansion  for  an  offset  TEMgg-mode  contains  thus  TEMf 
modes  only.  The  constants  are 


Am° 


2  2  1 

-e  /w 0  2  T  jo y  2  2 

e  (  2  m  )  /  dx{H.J  w  )exp^  "2(x  +  eX)/wo^>' 

W0lr2  m!  —  m  w° 


(11 


This  can  be  written  as 
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Amo  =  exp(-e2/2wo)(-e/wo)m/^' 


(1‘ 


Now  the  injected  electric  field  can  be  written  as  a  single  sum 


En  =  I  A_  E_ 


m=o 


mo  mo 


(1! 


or 


E0=/2nP  J— t 


e 

2W; 


2  2 
X  +j 


-  ^ 


W 


TTW 


e  exp[-2R“(x2+y2)]  E  [ (m! ) ” 1 ( -  ~ )m( 2 ) 

wo 


m=o 


u  r// 2xx-i* 

*  Hj(— )]x 


(1 


We  plot  equation  (17)  in  the  x  plane,  where  in  addition  to  z  =  0,  y 


i s  set  to  zero  as  wel  1 : 


2  2 
e  -x 


fo=/C2nP]  /  §-  e  e  I 

wn  tt  m=o  u 


(18) 


An  example  of  this  expansion  is  given  in  Figure  2  through  Figure 

.  2 

10.  The  dimensionless  amplitude  of  the  injected  beam  Eo=Eo(wo/8nP) 

vs.  (x/wo),  at  y=0=z  is  given  in  Figure  2.  The  beam  is  offset  by  an 

amount  e  which  is  set  equal  to  its  waist  size  wo.  The  result  of 

summation  for  the  first  seven  terms  is  given  in  Figure  10.  Figures  3  to 

9  give  an  impression  of  the  rate  of  convergence  for  this  choice  e  = 

w o •  This  is  an  enormous  misalignment.  Introducing  the  offset  parameter 

d  =  Z~  »  we  may  now  evaluate  the  magnitude  of  the  H-G's  created  by  an 
wo 

offset  TEM0g-mode  as 


A(Ci,m) 


=  -  (-d)' 

/[  m !  ] 


2 

-d  /2 


(19) 


II.  F.  2.  2)  Ti 1  ted  Source  Beam 

The  external  beam  crosses  the  resonator  at  its  waist  location  under 
an  angle  a.  This  is  shown  in  Figure  11. 

Again,  it  is  assumed  that  the  injected  beam  is  identical  to  the 
TEMgg-mode  of  the  resonator  except  that  it  is  rotated  at  the  waist  by  an 
angle  a.  The  electric  field  can  be  written  as 


E o  =  </C2nP]  (u '  00 


-jk(z  cos  a  +  x  sin  a) 


}(-z  sin  a  +  x  COS  a) 


,  2  ,z  _  r  ( -zsi  na+xcosct)  +y 

u  oo  =  {“I - 2 — F">  exP  l-  - 

w q  ir[l  +  (zcosa+xsina)  /z0]  Wg 


i  k  —  _ 

exp{-  [(-z  sin  a  +  x  cos  a)  +  y  ]}  e 


2  2 


j<f> 


The  vector  E0  can  be  written  as  follows. 


En  =  E  x  +  E  z 
u  xo  zo 


(21) 


where  Ex0  and  Ezo  are  the  x  and  z  components  of  the  vector  E0 
respectively.  Since  the  resonator's  eigenmodes  have  electric  fields  in 
the  x-direction  only,  see  Eq.  (4),  the  injected  electric  field's  x 
component  only  is  able  to  excite  the  resonator  modes.  The  y  component 
would  go  through  the  resonator  freely  as  there  is  no  conceivable 
coupling  of  an  orthogonal  vector  component  in  vacuum.  We  therefore 
expand  the  x  component  (i.e.  Exg)  of  the  field  E0  only: 


E 


xo 


E 

m=o 


Z 

n=o 


E 

mn 


(22) 


where 


mn 


/  /  E  E  dxdy 
‘  J  mn  xo  J 

-OO  -  00 


/  /  E  E  dxdy 
;  1  mn  mn  J 

-00  -oo 


Now  we  calculate  Amn  at  z  =  0  as  before. 


A 


mn 


f  f  dxdy{C  ip  ( - ) 

J  J  J  mnvnr Wn  1 


*  (-^) 
r  n  '  \a1  _  ’ 


rn' w 


}{(-2 - 

W07T[l+( 


xs i  na^  z-j" 
zo 


) 


2  2  2 

x  cos  a+y 
- T - 

w0 

-  e  cosa  e 


(-jkx  sin  a) 


(23) 


As  before,  it  is  possible  to  separate  the  double  integral 


_(X-CQ.S  q  +  jkxsina) 


mn  mn 


C  cos  a  .  /  d x {tp  — 

mn  j  wf 


zo 


2  2 


•  Jdy(*o(w„ 


(^)T» 


r 

TTWg 


(24) 


The  second  integral  is  zero  if  n  ^  0  or  it  is  unity  otherwise.  A^  may 
be  written  as 


mn 


V2x 


2  2 

xsi  na,  2t1/2  ,  x  cos  a 


Cmocosa  /  dx {4>m(— ^-)  [l+(^2)  ]'  exp(  ^ 


rnr  w 


n  *  0 

jkxsina)}  n  =  0 


w0  1  2  2  J’kw0 

Am0  =  Cm0  cos  a  ( - )  /  dC{H|T)(c)exp[-2  (1  +  cos  a)£  -(“Tj-sina)?] 

/2  _oo 


22  2 

Wg  £  sin  a  -  1  'z 

(1  +  - 2 - )  J 

2Zn 


(27) 


Amo  may  be  approximated  by 


V,  -  Cmo  cos  »  <77- >  ■’  d«  {Hm(5> 

-  00 

2  2 

exp[ - ( 1  +  cos  a)c  /2  -  jkwQsin  a  e//2]}. 


Let  us  calculate  the  integral  and  call  it  I. 


I 


/  dc{Hm(c)e“[a  C  +jSc]} 


where 


2  1+cos  a 
=  — 2 - 


kw0si na 

b  =  - 

✓2 


By  using  Euler's  formula,  "I"  can  be  written  as 


oo  _2  2  a,  .2  2 

I  =  /  d?[H  (e)e’a  c  cos(bc)]  -  j  /  d?[H U)e~a  c  sin  (be)] 


From  Ref.  [18],  we  have 


_2 

b 

I  .  4  e  ^  H  [ - ^ - 

3  ”  2a'2<77  - 

a 


l) 


■/2 


m 

( 1  -  ^ 
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Now,  we  can  substitute  for  a  and  b. 
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2n  ^  k  Wj  sin  a  kw0 


I  =  /( - iV)exp  [- 

1+cos  a 


7 — 3  •  [-, - -r-^1  [-  S1n  V3~ 

4(l+cos  a)  /2(1+C0S  a)  1+COS  a 


The  real  part  of  (i.e.  Re[Amg])  and  the  imaginary  part  of  A^g  ( 
Im[A  n])  are  written  as 


ReCAmol  Rjn 


for  even  in 


(33) 


ImCAm0]  =-jPm  for  odd 


where 


2  2  2 

; - k  w  0  sin  a 

P  =  jmwn  C  n  cos  a  / —  exp  [-  - 2 — ] 

™  0  m°  l*co77  ijl+coso) 


kwn/l+cos  a  ■ 

,,r  u  irSin  a  -iin 

HmC- - L-> - 

m  /2 (1+COS  a)  /1+COS  a 


where  Cn 


1  ,.m  . 

w  g  tt2  m! 


Now  we  substitute  for  Cm0  in  lm. 


P  =  j  cos  a  /- 


(1+cos  a)2  m! 


exp  [  - 


2  2  2 
k  w0  sin  a 

2  ^ 

4(1+C0S  a) 


nr  sin  a  ill 
rjl-r - 


/2(l+cos  a)  /1+cos  a 


In  order  to  plot  these  equations  and  make  a  numerical  example,  we 
would  like  to  inspect  the  function  vs  x  only  at  y  =  z  =  0  (waist).  The 
injected  beam  can  be  written  as 


Eo  =  x(  l  Pu)  +  z  E,„ 


at  z=o 


(36) 


Figure  12  shows  the  real  part  of  the  injected  field,  i.e. 


u00  =  u‘oo(x>  0»  0)  cos  (kxsina)  cosa  . 


A  tilt  given  by 


?2  =  a/tan"  (w0/L)  =  1 


is  chosen  is  the  tilt  angle  measured  in  terms  of  the  beam  spread.) 
For  a  waist  size  of  Wg  =  0.002m  =  2mm,  and  a  wave  number  k  =  (2ir/633 
xl0~9)m“l  =  9.93x10^,  this  corresponds  to  a  tilt  angle  of  a=1.0xl0-^  rad 
=  21  arc  second.  Fig.'s  13  and  14  show  two  components. 

In  figure  15,  the  sum  of  the  first  two  terms  is  shown.  Comparing 
it  to  fig.  12,  the  quick  convergence  is  evident.  Figure  16,  with  four 
terms,  is  indistinguishable  from  fig.  12. 

II.  F.  2.  3)  Mismatched  Source  Beam 

The  injected  TEMgg-beam  shall  have  a  waist  W.  The  resonator's 
eigenmodes  shall  be  Hermite-Gaussian  modes  with  waist  size  wo.  At  z  = 
0,  we  expand  the  injected  beam  in  terms  of  these  Hermi te-Gaussians ,  (see 
Figure  17). 

The  injected  electric  field  is 

r _ -Ji[xV] 

to  ■  /V-  Ae"  e  *  .  eJ*  e'Jk2  ‘x  .  (39) 


As  before,  we  expand  EQ  in  the  plane  z  =  0. 


E0=  /2nP-  2/irW0  exp[-(x  +  y  )/W0]. 


Finally  the  injected  electric  field  can  be  written  as  the  following 
summation. 


2  2 
x  +y 


x  +y 


"° 


w0  +W0 

n*°  ,2,. .  %  ?!  2"*"  (M0%  $  1X 


A  numerical  example  is  given  in  Figures  18-24.  In  mismatching  problems 
as  given  above,  we  may  define  a  mismatch  parameter 


!  Wq  +  W Q 

In  Figure  18,  the  injected  beam  amplitude  is  shown,  in  arb.  units.  It 
has  a  mismatch  £3  =  0.77  (Wq  =  2wq). 

Figuresl9  to  22  show  Uoo,  U02  ,U20  and  Irrespectively. 

In  figure  23,  the  summation  of  the  first  2x2  terms  is  shown: 

Uoo  Aoo  *  U02  A20  +  U20  A 20  +^,12.  • 


Ei'aen  modes 


Figure  23  shows  a  kurtosis.  In  Figure  24,  the  summation  of  the  first  8 
x  8  terms  is  shown. 
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m=°  2 

hi  ,  ,  .  . 
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mn 
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mn 


The  peak  of  this  partial  sum  is  within  99%  of  the  peak  of  the  original 
curve. 

II.  F.  3.  Application  of  the  results  of  section  II.  F.  to  ringlasers 

The  mi sadjustment  of  an  external  stigmatic  gaussian  beam  to  a  ring 
gives  rise  to  Hermite-Gaussian  eigenmodes  in  addition  to  excitation  of 
the  ring  in  the  proper  Gaussian  mode.  An  added  complication  is  the 
presence  of  astigmatism  of  the  ring  modes.  The  latter  can  be  expressed 
by  the  presence  of  an  elliptic  ring  waist  with  the  two  waist  sizes  wQX 
and  woy. 

A  misalignment,  in  two  dimensions,  is  then  expressed  by  the  offsets 
ex  and  ey,  and  the  tilt  angles  ax  and  ay,  see  figure  25.  A  slight 
generalization  of  the  previous  section  gives  then  the  misalignment 
parameters 


e  /w 
x  ox 


{ q  /w  )  a 
^ox  ox  x 


Cy  "  ey/woy  +  ^oy^oy^y  ’ 

where  qox,  qoy  are  the  complex  curvatures  of  the  ring  modes  at  the 
waists.  This  simple  set  of  equations  is  proper,  since  the  propagation 
characteristics  of  the  beam  in  the  two  orthogonal  planes  are  not 
coupled. 


■”W'- 


Furthermore,  a  mismatch,  see  figure  26,  is  given  by  the  mismatch 
parameters 

n  ,w°  '  wox,  V2  n  _  '  wov//2 

*  (“2 - 7~ )  »  °v  -  (“2 - T~ I 

x  "  -  -  y  W«  +  w_ 


W  n  +  W. 


Note  that  for  perfect  adjustment,  all  parameters  are  zero. 

The  power  ratio  of  a  H.-G.,  Pmn/Poo  to  the  power  of  the  Gaussian  in  the 
ring  can  then  be  expressed  as 


^mn^00  =  *  ^n^y’^y^’ 


F„<V0x>  '  <2mm!)-1|0x|2m  I H Jlj < l-0*)«x  . 


F „(5y,0y)  *  (2"n!)'‘|0y|2n|Hn[j(l-0y)Ey  //20y|2  , 


judiciously  exploiting  the  previous  sections.  These  equations  are 
useful  to  estimate  the  amount  of  pulling  due  to  offset,  misalignment, 
and  mismatch. 

The  power  transfer  from  a  misadjusted  Gaussian  to  the  ring’s 
Gaussian  is,  of  course,  of  interest  as  well:  For  pure  mismatch,  with 
coinciding  waists,  it  is  simply 

p o o/P|_  =  1  -  04  (0  ->■  0X  or  0y,  ?i  =  Power  of  injecting  laser). 

For  pure  offset, 

P  o  o/P{_  =  exp(-£2/wo)  =  exp[-(Re5)2] 

The  figure  27  shows  a  plot  of  the  latter  two  equations.  They  show  that 
alignment  is  to  be  treated  more  carefully  than  matching. 


/og 


II.  F.  4  Experimental  test  of  misalignment  and  mismatch. 

The  previous  notions  were  put  to  the  test  by  applying  them  to 

experimental  results  obtained  by  the  Seiler  group  (ref.  12). 

In  a  ring  injected  by  a  presumably  stigmatic  Gaussian  beam  with  no 
matching  circuit,  the  scanning  analysis  produced  figure  28.  The  two 

largest  spikes  were  identified  as  q,0,Q  and  q+1,0,0  modes,  the  rest  were 
labelled  "higher-order  modes".  An  analysis  of  the  frequencies  involved 
verified  that  they  are  all  H.-G.'s,  with  indices  as  given  in  the  figure. 
15  such  modes  were  identified. 

The  fact  that  one  of  the  most  prominent  H.-G.  was  qOl  indicated 

serious  mi  sal ignment,  as  mismatch  produces  even-indexed  modes  only. 

The  fact  that  the  qlO  mode  is  not  visible  indicates  alignment  in 
the  plane  of  the  ring. 

The  fact  that  the  q02  mode  and  the  q04  mode  were  almost  equally 

strong  indicates  misalignment  in  the  y-direction  in  general  and  no  major 
mi  smatch. 

On  the  other  hand,  the  presence  of  the  q20  and  the  q22  modes  allows 
an  estimate  of  the  mismatch  in  the  x-direction. 

Given  the  optical  circuit,  the  mismatch  parameters  are 
0X  =  -0.53  ±  j0.66  ,  Oy  -  -0.40  ±  j0.76 

with  |0X |  =  0.84  ,  |0y |  *  0.87 

The  misalignment  parameter  in  the  x-direction  is  negligible  (good 
alignment  in  the  plane  of  the  ring),  but  in  the  y-direction  it  is 
=  0.255  -  j 0.1 53 

as  obtained  from  a  least-squares  fit  of  the  H.-G.  powers  of  lines  qOl 
through  q012  to  the  equation  for  misalignment  (The  fit  was  good  to 
2.4%). 


The  above  result  says  that  the  offset  in  the  y-direction  was 
ey=180um,  and  the  tilt  in  the  y-direction  was  ay  =  -49  arc  second. 

The  same  paper  shows  a  "partly  matched"  beam.  Fig.  29,  where  a 
spherical  lens  was  placed  in  a  compromise  position.  The  pattern  shows, 
besides  the  q02,  q20,  and  qOl  modes,  a  newly  emerging  qlO  mode 
indicating  that  now  there  is  also  misalignment  in  the  x-direction. 
the  parameters  are 

0X  =  0.328  +  jO,  0y  =  0  +  j0.303,  ||x|=0.16,  |cy|=0.19 
reflecting  a  bigger  effort  in  alignment  as  well  as  match. 

Finally,  the  "matched"  (actually  aligned)  cavity  shows  the 
fundamental  Gaussians  only,  figure  30. 

For  "tutorial"  purposes,  figure  31  shows  the  cross  sections  of  an  offset 
and  mismatched  beam  and  its  expansion  into  all  Hermite-Gaussian 
eigenmodes  of  the  astigmatic  resonator.  The  contours  are  all  the  same 
isophote.  Also,  the  aligned  laser  is  shown  to  produce  only  even-indexed 

H. -G.'s.  Note  that  the  width  Wg  of  the  circular  incoming  beam  is  made 

equal  to  the  horizontal  width  wQX  of  the  resonator. 

These  results  suggest  an  interesting  way  to  keep  a  large  ring 
aligned  via  servos: 

I.  Pick  off  the  qlO  mode  and  control  with  it  the  horizontal  injection. 

2.  Pick  off  the  qOl  mode  and  control  with  it  the  vertical  injection. 

3.  Pick  off  the  q02  mode  and  control  with  it  the  vertical  complex 
curvature  of  the  beam  (y-zoom). 

4.  Pick  off  the  q20  mode  and  control  with  it  the  horizontal  complex 
curvature  of  the  beam  (x-zoom). 


III.  A.  Scanning  of  an  astigmatic  beam,  and  beam  evolution: 

New  method  to  find  beam  and  evolution  parameters 

This  section  contains  the  first  experimental  subject  we  were 
working  on.  It  was  motivated  by  the  desire  to  establish  experimentally 
the  shape  and  size  of  an  astigmatic  beam  spot.  Scanning  by  a  razor 
blade  (“knife-edge")  was  recognized  to  deliver  enough  information  to 
find  the  axes  of  an  elliptical  spot,  and  its  inclination.  The  analysis 
proved  to  be  surprisingly  simple,  and  its  results  easily  applicable. 
Data  evaluation  via  non-linear  least  squares  fitting  is  applied  to  find 
spot  parameters  as  well  as  beam  evolution.  The  programs  are  available 
upon  request. 

In  the  following,  a  reprint  of  a  paper  is  given  (ref.  19)  which 
summarizes  this  work. 


Knife-edge  scanning  of  an  astigmatic  Gaussian  beam 


Hans  R.  Bilger  and  Taufiq  Habib 


The  relations  for  position,  spot  size,  and  inclination  of  the  major  axis  of  an  elliptical  Gaussian  beam  to  knife- 
edge  scanning  data  are  derived.  A  knife-edge  whose  scanning  direction  is  adjustable  to  any  angle  has  been 
employed  to  scan  across  a  beam  in  at  least  three  directions.  Nonlinear  least-squares  fit  programs  have  been 
developed  to  check  whether  a  beam  is  Gaussian,  and  to  evaluate  the  parameters,  with  errors,  of  such  an  ellip¬ 
tic  spot.  The  evolution  of  an  astigmatic  beam  in  the  tangential  and  sagittal  plane  is  measured. 


I.  Introduction 

Scanning  of  a  Gaussian  beam  with  circular  cross 
section  has  been  dealt  with  in  several  papers.1  How¬ 
ever,  beams  with  different  cross  sections  or  power 
density  distributions  appear  often  in  laser  systems  ei¬ 
ther  as  aberrations,  e.g.,  through  admixture  of  Her- 
mite-Gaussian  beams,  or  as  an  essential  feature,  as  is  the 
case  in  ring  lasers.  The  extension  reported  here  has 
several  benefits  in  the  laboratory: 

(1)  Quantitative  determination  of  the  degree  to 
which  a  given  spot  deviates  from  a  circular  cross  sec¬ 
tion. 

(2)  Determination  of  the  parameters  of  an  elliptical 
spot. 

(3)  Analytical  (least-squares)  evaluation  of  the  spot, 
including  an  estimate  of  the  errors  of  the  evaluated 
parameters. 

II.  Scanning  of  a  Gaussian  Beam  with  Elliptical 
Cross  Section 

The  spot  size  w  shall  be  defined  by  the  following 
distribution  of  the  electric  field  amplitude  E: 

£(x,y)  *  Eo  exp[-(x2  +  y2)/w2\  (1) 

for  a  circular  beam  spot  centered  at  x  =  y  =  0;  the  beam 
propagates  in  the  direction  of  the  positive  z  axis. 

Extending  this  to  an  elliptical  cross  section  with  ar¬ 
bitrary  orientation  (see  Fig.  1)  and  using  power  density 
S  (proportional  to  intensity  /),  we  have 

2  P 

S(x'y)  ■ - —  exp(-2x'2/uij)  exp(-2y'2/u/|),  (2) 

TWaWt, 
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where  wa,wb  are  the  spot  sizes  along  the  major  and 
minor  axes,  respectively,  P0  is  the  total  power  in  the 
beam,  2(PoJirwawb)  is  the  power  density  in  the  beam 
center,  and 

x  -  xo  *  *'  cosao  —  y'  sinao. 

y  ~yo  *  x’  sinaro  +  >'  cosao-  (3) 

Equations  (3)  include  rotation  by  a  tilt  angle  a0  (mea¬ 
sured  from  the  positive  x  axis  in  the  direction  of  the 
positive  y  axis)  and  translation  of  the  beam  center  to 
x0,yo-  A  centered  circular  Gaussian  beam  [Eq.  (l)j  is 
included  as  a  special  case  with  x0  =  yo  =  0  and  wa  *  wb 
-  w. 

When  a  straightedge  is  placed  at  x  =  xs  parallel  to  the 
y  axis,  obscuring  the  half-plane  x  $  xs,  the  power 
transmitted  past  the  edge2  is  given  by 

/’(X.)  -  £7  lx:::  S(x,y)dvj  dx.  (4) 

This  double  integral  can  be  expressed  as  a  comple¬ 
mentary  error  function  of  xs,3  namely, 

P(x.)  -  (Po/2)  erfc(u),  u  -  (v/2)(x.  -  x,)/u><ao),  (5) 

with  erfc(u)  =  (2 /•>/*■)/"  exp(-t2)dt,  or  as4 

P(x,)/Po  «  Q(v),  v  »  [2/ui(ao)l(x,  -  x,)  *  (v^2)u,  (6) 

with 

Q(v)  -  [lV(2ir)|  j *"  exp(-  tV2)dt, 
and 

wHcio)  “  w\  co«Jao  +  u’l  sin2ao, 

x,(a0)  *  xo  cosa0  +  yosin«o  +  x*,  (7) 

(x,0  is  the  offset  of  the  translator  micrometer). 

The  function  P(x,)/P0  is  drawn  in  Fig.  2  vs  scanner 
position  x,  for  two  widths,  u;i  and  ic2;  the  center  of  the 
beam  is  assumed  to  be  at  x,  *  0.  Equation  (7)  contains 
the  three  parameters  of  interest  w„,  wb,  and  «o-  The 
apparent  width  w(ao)  is  plotted  in  a  polar  plot  in  Fig. 
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Fig.  1.  Elliptical  beam  spot  centered  at  i<>,yo  with  axes  wQ  and  wt, 
tilted  against  the  positive  x  axis  (horizontal)  by  oq.  The  beam  goes 
into  the  paper  plane  (positive  z  axis).  The  scanning  edge  is  placed 
at  xa. 


Fig.  3.  Polar  plot  of  widths  vs  scanning  angle  a  for  three  different 
elliptical  beam  spots  (solid  curves).  The  spot  for  wa/u'b  =  2  is  drawn 
as  a  dashed  ellipse. 


* 


Fig.  2.  Relative  power  P(i,  )/Po  vs  scanner  position  x,  for  two  spot 
sizes  w  i  ( »  1.5  units  on  the  abscissa)  and  iuj  (»  3  units  on  the  abscissa). 
The  shapes  of  these  curves  are  identical,  except  for  a  difference  in  w, 
for  scanning  of  an  arbitrarily  placed  elliptical  beam.  The  beam 
centers  are  assumed  to  be  at  x,  >0. 

3  for  three  different  ratios  wa/wb,  together  with  the 
actual  spot  shape  (ellipse,  dashed  line)  for  wa/wb  *  2. 
The  measured  widths  w(ao)  agree  with  the  widths  of  the 
ellipse  at  a  *  a0  and  a  *  a0  ±  180°  where  u;(ao)  *  wa, 
and  at  a  ■  a0  ±  90°  where  w(ato)  =  w 
If  the  scanner  is  now  rotated  around  the  z  axis  by  an 
angle  a  against  the  positive  x  axis,  the  equation  for  the 
resulting  width  is  slightly  generalized  to 

u/2(a)  ■  uil  cos *(a0  -  a)  +  u>\  sin2(ao  -  a), 

(8) 

x,(a)  »  xocostao  -  or)  +  yosin(ao  -  a)  +  x,0. 

For  a  circular  beam  with  w„  m  wt,  *  w,  the  result  is  again 
w(a)  m  const  *  w,  as  it  should  be. 


Three  methods  have  been  devised  to  evaluate 

U)(a): 

(a)  Fractional  power  method:  Noting  that  Q(v+  = 
1)  =  0.841  =  P(x*)/P0 ,  and  Q(u-  =  -1)  =  0.159  = 
P{xJ)/Pq  (see  Fig.  2),  we  set  the  edge  such  that  the 
relative  power  equals  these  fractions.  The  positions  x  f 
and  x~  yield  the  width  w(a)  through 

-  v-  »  2  «  [2/u;(a)](x*  -  x.)  -  [2/u)(a)|(x,"  -  x,), 

or 

w(a)  =  x,»  -  x;.  (9) 

In  principle,  any  pair  o+o_  can  be  chosen  to  evaluate  w, 
but  with  a  given  absolute  error  in  P(v ),  it  can  be  shown 
that  the  pair  u±  -  ± 1  produces  the  width  with  mini¬ 
mum  relative  error.  The  center  of  the  beam,  xs ,  can  be 
obtained  by  setting  the  scanner  such  that  the  trans¬ 
mitted  power  is  halved,  see  Fig.  2. 

(b)  Graphical  method:  This  method  consists  of 
plotting  the  relative  power  P(x,)/P0  vs  x,  on  error 
function  paper.5  A  straight  line  on  this  paper  indicates 
an  error  function,  i.e.,  that  the  beam  is  indeed  Gaussian. 
The  points  x*  and  x,“  can  then  be  used  as  above  to  find 
w.  The  center  of  the  beam  is  again  given  by  P(x,  )/P0 
=  1/2.  This  method  has  the  advantage  over  method  (a) 
in  that  it  makes  use  of  all  the  measured  scanning  posi¬ 
tions,  that  the  latter  do  not  have  to  be  specifically  cho¬ 
sen,  and  that  the  graph  allows  a  check  whether  the  beam 
is  Gaussian.  It  also  provides  some  estimate  of  mea¬ 
surement  errors. 

(c)  Least-squares  method:  This  third  method  con¬ 
sists  of  fitting  the  Q- function  into  the  measured  power 
vs  x,  (see  Fig.  2),  by  adjusting  the  three  parameters  Po, 
w,  and  x,.  This  method  does  not  require  the  data  to  be 
normalized  with  the  power  Po  before  analysis.  It  can 
furthermore  accommodate  a  fluctuating  total  power  and 
it  reduces  errors  introduced  by  fluctuating  power.  The 
accuracy  of  the  computer  program  is  not  limited  to  the 
usual  ~1%  of  graphical  methods.  Any  observer  bias  is 
eliminated.  Finally  the  program  also  provides  analytic 
estimates  of  the  errors  of  the  parameters  by  calculating 
the  variance-covariance  matrix.6 
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Fig.  4.  Test  of  a  spot  at  z  *  90  cm  off  the  laser  without  intervening 
optics  at  5  angles,  0",  45®,  90°,  135°,  and  180®.  The  fit  ui(a)  “  const 
results  in  w  »  (653  ±  4)  Mm.  The  sensitivity  of  the  method  to  detect 
ellipticity  is  demonstrated  by  the  dashed  curve  where  tentatively 
wjwa  «  0.9  has  been  set. 


Given  the  measured  points  P(xs),  usually  between 
five  and  twenty  points  covering  the  range  from  0.9  Po 
to  0.1  Po,  the  program  fits  the  function 

P(z)  =  alQ(2(x  -  a3)/02l  »6a) 

into  the  points  where  at  =  P0,  <12  =  w(a),  03  =  x,.  An 
approximation  for  Q  is  chosen4  which  has  a  maximum 
error  of  ±  1  X  10"5  over  all  arguments  v: 

Q(v)  ~  +  sgn(o)|[l/v,'(2*-)| 

X  (a  +  6r  +  ct2)t  exp(-u2/2)  -  1/21, 

t  -  1/(1  +p|H>. 

sgn(u  <  0)  * -1,  sgn(u  =  0)  =  0,  sgn(o  >  0)  »  +1, 
p  «  0.33267,  a  =•  0.4361836, 
b  •  -0.1201676.  c  -  0.9372980. 

The  program  is  iterative.  It  makes  use  of  initial  esti¬ 
mates  for  o i,  02, 03  and  refines  these  estimates  until  the 
sum  of  the  residuals  squared  is  satisfactorily  close  to  the 
minimum.  For  a  set  of  twenty  points,  the  program 
takes  ~800-msec  CPU  time  on  the  VAX  11/750  for  four 
iterations,  which  usually  leads  to  convergence.7 

III.  Parameters  of  a  Gaussian  Beam  with  Elliptical 
Cross  Section 

One  scan  at  an  angle  a  gives  w(a )  =  02-  We  need  at 
least  three  scans  to  calculate  wa,Wb,  and  ocq-  In  prac¬ 
tice,  scans  at  more  than  three  angles  are  made.  A  sep¬ 
arate  least-squares  fit  program  then  determines  the 
three  parameters  above  through  Eq.  (8)  or  through 

w(a)  m  w,  v'll  -  «2sin2(ao  -  a)|,  (8a) 

with  t  =  numerical  eccentricity  *  Vd  “  w\/w\). 

IV.  Establishment  of  an  astigmatic  Gaussian  beam 

A  He-Ne  laser  (Oriel  model  6697)  was  used  to  pro¬ 
duce  a  well-behaved  circular  Gaussian  beam.  Its  cir¬ 
cularity  was  checked  by  scanning  at  several  angles,  see 
Fig.  4.  A  silicon  detector  (Optics  Technology  model 
610)  with  a  narrowband  optical  interference  filter  was 


Fig.  5.  Goniometer  to  adjust  the  obliquity  angle  of  the  lens.  The 
reference  position  is  given  by  retroreflection  which  originally  aligns 

the  principal  axis  of  the  lens  with  the  beam  axis  =  optical  axis. 

used  for  detection  with  background  light.  The  beam 
evolution  was  measured  within  the  first  1.5  m  off  the 
laser  front  end.  The  spot  size  vs  z  was  fitted  by  the 
equation 

w(z)  -  w0\  [1  +  (x  -  zw)2/zl\,  (11) 

with  Wo  =  waist  size,  located  at  z  =  zw,  and  z0  =  Ray¬ 
leigh  length  =  ttwI/X.  The  waist  location  zw  was  found 
near  the  output  mirror  of  the  laser.  The  waist  size  was 
wq  -  (316  ±  5)  nm. 

A  plano-convex  lens  with  focal  length  fo  =  20  cm  was 
then  placed  on  a  rudimentary  goniometer  (Fig.  5),  which 
enabled  us  to  rotate  the  lens  around  a  vertical  axis  by 
an  angle  X*  and  also  around  a  horizontal  axis  by  an 
angle  X„.  The  angle  of  obliquity8  <t>  between  the  beam 
and  the  principal  axis  of  the  lens  then  becomes 

<)>  »  arccos(co8>*  cosy,,).  (12) 

This  angle  lies  in  the  tangential  plane  whose  tilt  against 
the  horizontal  plane  (x  -z  plane)  is  given  by 

a0  =  arcsin(sin7„/sim#).  (13) 

The  tangential  plane  therefore  contains  one  of  the  axes 
of  the  ellipse. 

The  lens  was  placed  at  2/0  =  40  cm  from  the  beam 
waist,  with  the  two  angles  -  yv  «=  30°,  which  results 
in  the  two  focal  lengths9:  /unf*nu»i  -  9-4  cm  and 
/sagittal  =  16.8  cm.  The  beam,  after  traversing  the  lens, 
is  expected  to  have  a  shape  as  given  in  Fig.  6,  with  the 
circle  of  least  confusion8  9  at  z  *  15.5  cm;  before  and 
after  this  point,  the  ellipse  rotates  by  -90°  (a0  *  136.5° 
to  ao  =  46.5°).  Since  the  input  to  the  lens  is  a  circular 
Gaussian  beam,  there  will  be  a  tangential  waist  and  a 
sagittal  waist  placed  approximately  symmetrical  to  z<iC. 
There  are  no  focal  lines  in  this  case. 

The  propagation  of  the  beam  in  the  tangential  plane 
is  independent  of  that  in  the  sagittal  plane.  Both  are 
governed  by  Eq.  (11),  with  two  different  sets  of  pa¬ 
rameters:  Wot,  zwt,zot.  and  wo,,  zw,,zo,. 

V.  Experimental  Evaluation  of  the  Astigmatic  Beam 

The  beam  was  probed  at  distances  12-60  cm  after  the 
lens.  This  range  contains  the  interesting  features  (see 
Fig.  6):  primary  focus  =  waist  in  the  tangential  plane, 
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Fig.  6.  Evolution  of  beam  along  the  2  axis.  The  circular  Gaussian 
beam  enters  the  lens  at  r  =  0.  The  spots  are  shown  as  they  appear 
in  the  x-y  plane  enlarged  by  x  100  relative  to  the  z  scale.  Immediately 
to  the  left  and  right  of  the  circle  of  least  confusion,  the  tangential  and 
sagittal  waists  are  shown,  respectively. 


Fig.  7.  Scan  of  astigmatic  beam  at  2  «  60  cm  from  the  lens  at  or  =  45°. 
The  least-squares  fitted  parameters  are:  total  power  P0  =  (1.498  ± 
0.005)  mW,  center  at  x,  =  (6.309  ±  0.002)  mm,  spot  size  u;(450)  =  (783 
±  6)  «im. 

circle  of  least  confusion,  and  secondary  focus  =  waist 
in  the  sagittal  plane. 

Figure  7  shows  a  typical  measurement  of  P(x)  at  a 
distance  z  -  60  cm  from  the  lens,  at  a  scanning  angle  a 
-  45°  with  respect  to  the  positive  x  axis.  The  least- 
squares  fit  of  Eq.  (6a)  to  the  fourteen  measured  points 
gives  P0  =  (1.498  ±  0.005)  mW,  u;(450)  =  (783  ±  6)  /an, 
x,  =  (6  309  dk  0.002)  mm.  None  of  the  individual  re¬ 
siduals  exceeded  8  /tW;  a  translator  with  10-^m  reso¬ 
lution  (smallest  division)  was  used. 

After  at  least  four  widths  are  evaluated  for  each  spot, 
Eq.  (8a)  is  fitted  into  the  data  w(a).  Figure  8  shows  a 
typical  fit  at  z  *  35  cm  with  the  result  u/a  =  (801  ±  3) 
Mm,  c  *  0.920  ±  0.003,  a0  =  (46.4  ±  0.4)°,  which  estab¬ 
lishes  the  tilt  and  size  of  the  ellipse  at  35  cm  together 
with  the  errors.  The  maximum  deviation  is  7.5  fim;  the 
average  deviation  is  4  nm. 

Finally,  the  evolution  of  the  tangential  spot  sizes 
w,(z)  and  the  sagittal  spot  sizes  w,(z)  is  obtained  (see 


Fig.  8.  Polar  plot  of  spot  sizes  w(a ),  taken  at  z  *  35  cm.  The  fit  gives 
u>q  =  (801  ±  3)  nm,  t  =  0.920  ±  0.003,  a0  -  (46.4  ±  0.4)°.  The  rms 
deviation  of  the  widths  in  this  plot  from  the  best-fitted  solid  curve 
is  4.5  11m.  The  resulting  elliptical  spot  size  is  drawn  as  a  dashed 
curve. 


Fig.  9.  Evolution  of  astigmatic  beam  vs  z  (see  also  Fig.  6).  The  waist 
sizes  are  ieo»  «  (58.0  ±  0.2)  |im  and  wo.  »  ( 102.0  ±1.2)  urn.  They  can 
be  located  with  an  accuracy  of  about  ±2  mm. 


Fig.  9)  by  least-square  fitting  Eq.  (11)  into  the  previ¬ 
ously  obtained  results.  The  best-fitted  parameters 
are 

tangential  waist  tMor  *  (58.0  ±0.2)  Mm  atZu,,  =  (118.5 
±  0.9)  mm, 

sagittal  waist  u>o$  =  (102.0  ±  1.2)  Mm  at  zws  ~  (201.3 
±  2.5)  mm. 

The  maximum  deviation  of  any  point  was  8.5  nm  for 
the  tangential  spot  sizes  and  14.6  Mm  for  the  sagittal 
spot  sizes. 

VI.  Discussion  of  Results 

The  methods  used  here  to  scan  an  elliptical  Gaussian 
beam  with  a  straightedge  (razor  blade)  are  based  on  the 
result  that  the  power  spilling  over  the  blade  has  the 
same  dependence  on  the  position  of  the  blade  as  is  the 
case  in  a  circular  Gaussian  beam,  except  that  the  eval¬ 
uated  width  is  a  simple  function  of  position  and  size  of 
the  ellipse.  At  least  three  scans  need  to  be  done  on  each 
spot,  with  different  angles  a.  If  the  position  of  the 
tangential  plane  («0)  is  known,  two  measurements 
would  suffice,  preferably  taken  at  ct  ■  a0  and  a  *  a0  ± 
90°. 
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Two  methods  were  mainly  used  to  obtain  a  width, 
namely,  the  fractional  power  method  and  the  least- 
squares  method.  The  latter  established  that  the  beam 
is  indeed  Gaussian,  within  the  errors.  However,  sub¬ 
sequent  measurements  were  usually  done  with  the 
fractional  power  method,  which  is  faster.  Both  meth¬ 
ods  yield  typical  errors  of  the  order  of  one  small  division 
on  the  scanner  (=  10  /im)  or  less.  The  data  suggest 
indeed  that  translators  with  a  resolution  of  1  would 
produce  yet  smaller  random  errors. 

To  test  the  sensitivity  of  the  methods  to  detect  de¬ 
viations  from  Gaussian  profiles,  a  Hermite-Gaussian 
1-0  mode  of  the  same  spot  size10  was  evaluated.  The 
calculated  result  of  scanning  along  the  x  axis  is 

(Pfx.J/PolrEM.,,)  *  (1/2)  erf c(u)  +  (u/v/ir)  exp(-u2), 

,  _  (14) 

u  =  (v^2)(x  -  x,)/w 

The  major  deviation  from  the  error  function  occurs  at 
um  =  ±1/V2  with  an  amount  of  y/(2ire)  =  24%.  This 
suggests  that  the  detection  of  admixtures  of  such  ei- 
genmodes  with  a  power  of  less  than,  say,  10%  of  the 
fundamental  Gaussian  has  to  make  use  of  other  meth¬ 
ods,  e.g.,  of  Fabry-Perot  scanning  in  the  frequency  do¬ 
main. 

It  may  be  noted  that  by  using  the  least-squares  pro¬ 
gram  above,  an  iteration  of  the  type  given  in  Ref.  11  can 
be  avoided:  The  program  finds  directly  one  waist  to 
satisfy  both  the  asymptotic  slope  l/z<>  =  A/ira'o  as  well 
as  the  minimum  woof  the  function  to  (2 )  =  u)o\/U  +  (2 
-  zw)2/zl\  with  the  criterion  of  minimizing  random  er¬ 
rors. 


This  paper  is  an  outcome  of  research  done  under 
AFOSR  grant  84-0058. 
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III.  8.  Vacuum  in  ring:  Residual  Fresnel  drag 
A  passive  ring  can  be  evacuated.  What  is  the  required  residual 
pressure?  We  investigated  two  effects  on  the  output,  here  we  report  on 
Fresnel  drag  (section  III.  C.  shows  calculations  on  quality  factor). 

A  portion  of  this  section  is  published  (ref.'s  20,  21). 


r*  < 


Vacuum  in  ring:  Residual  Fresnel  drag 
Light  drag  is  a  surpri singly  simple  means  to  introduce  an 
anisotropy  in  the  ring.  Experimental ly ,  some  drag  equations  have  been 
verified  to  the  level  of  0.1%.^ 

For  a  homogeneous  velocity  field  v  of  a  gas  parallel  to  the  light 
beam  over  a  distance  d,  the  beat  frequency  in  a  ringlaser  is  to  a  good 
approximation  given  by 

Af  =  (4/xL)  (n-1)  vd,  (1) 

where  n  =  index  in  refraction.  For  helium  this  approximation  is  good  to 
1%.  The  quantity  (n-1)  is  proportional  to  the  gas  pressure.  For  an 
ideal  gas 


(n-i)  =  (T0/T)  (P/Po)  (n-l)T  n 

1 > A  1 o  »Po 


(2) 


T  =  absolute  temperature  of  flowing  gas,  p  =  pressure  of  flowing  gas. 
To,  po  =  standard  conditions.  For  stationary  mass  flow,  the  flow  rate 
M  =  d(pV)/dt  =  q(vp)  =  kT  dN/dt  =  constant  (3) 

(q  =  tube  cross  section,  V  =  volume  of  gas,  N  =  number  of  gas  molecules, 
k  =  Boltzmann  constant).  It  follows  for  the  beat  frequency  due  to  light 
drag 


Af  =  (4/aL)  (d/q)  ( n-1 )T  . 

1 0>Po  »A 


(T0/Tp0)M 


(4) 


The  beat  frequency  is  therefore  determined  by  the  mass  flow  rate.  The 
velocity  vs  length  is  arbitrary  and  need  not  be  known;  the  drag  is 
independent  of  pressure  and  pressure  gradients  along  the  path,  as  long 
as  M  is  constant  and  the  flow  is  stationary. 


Vacuum  Requirement.  Inadvertent  drag  due  to  moving  residual  gas  has  to 
be  avoided.  From  the  foregoing  we  can  estimate  upper  limits  of  mass 
flow  sources  for  negligible  drag.  To  push  drag  effects  below  10"^  Hz 


generally  we  require  the  sources  (inadvertent  leaks,  outgassing,  etc.) 
to  have  flow  rates  below  10*^  Pa*m8/s,  which  is  equivalent  to  10"8  torr 
£/s. 

However,  typical  information  about  vacuum  is  given  through 

pressure,  which  in  turn  is  related  to  flow  via  the  pump  resistance  of 
the  ring. 

It  turns  out  that  with  mass  flow  rates  even  as  high  as  "100  micron 
cubic  foot  per  hour"  (10-^  Pam8/s)  we  are  already  close  to  the  high 

vacuum  regime.  In  this  regime,  the  pump  resistance  R  (s/m8)  of  a 

cylindrical  tube  of  cross  section  q  and  length  d,  molecular  weight  m  is 

R(s/m3)  =  ( 3 it /4 )  (m/2kT)1/2  d/q3/2  (5) 

(m=mass  of  one  molecule)  and  the  pressure  difference  along  such  a  tube 
is 

P2  -  Px  =  MR  .  (6) 

Assuming  a  worst  case  of  "inadvertent"  mass  flow  injection  of  M  =  10“^ 
Panr/s  farthest  away  from  the  pump,  d  =  7.6  m,  and  a  tube  with  1  cm 

diameter  as  above,  we  obtain 

P2  -  Px  =  P2  =  0.62xl0"4  Pa  (=4.7xl0-7  torr) 

whereby  the  pressure  P^  at  the  input  of  a  75  1/s  pump  would  be  1.3xl0-8 
Pa  (=  lxlO’10  torr),  i.e.,  negligible.  From  this  we  conclude  that 
reasonable  high  vacuum  conditions  have  to  prevail.  They  are  achievable 
in  a  laboratory  with  modest  cleaning  and  pumping. 

It  should  be  noted,  that  an  "active"  ringlaser  requires  by 
definition  material  in  the  light  path,  e.g.  the  HeNe  gas  for 


amplification.  Indeed  much  of  the  residual  problems  with  small 
commercial,  active,  ringlasers  are  directly  related  to  light  drag  from 
moving  gas  (Langmuir  flow)^.  The  gas  pressures  needed  to  sustain  the 
necessary  amplification  are  so  far  above  those  calculated  in  this 
section,  that  it  is  considered  impossible  to  render  negligible  any  false 
signals  (noise)  due  to  Fresnel  drag.  We  consider  this  true  even  if 
elaborate  measures  are  used  to  balance  mass  flow  etc,  as  it  is  done  in 
small  ring  laser  gyros. 


IZ 

III.  C.  Vacuum  in  ring:  Effect  on  quality  factor 

Power  loss  of  a  beam  due  to  scattering  by  a  gas  was  recognized  very 
early  as  a  source  of  imperfection,  but  it  was  clearly  considered 

negligible.  After  more  important  error  sources  had  been  eliminated, 

this  effect  got  more  into  the  forefront.  In  a  large  ringlaser,  with 
very  high  sensitivity  and  quality  factors  of  the  order  of  10^, 

scattering  becomes  a  substantial  effect.  Fortunately,  as  this  section 
demonstrates,  even  modest  pumping  will  render  negligible  its  effect  on 
the  quality  factor  of  even  a  large  resonator.  Rayleigh  scattering  off 
neutral  atoms  is  treated  as  the  only  scattering  mode.  Born  and  Wolf‘d 
served  as  a  basis  of  a  classic  (non-quantum)  calculation  in  this 


section. 


Vacuum  in  a  ring:  Effect  on  quality  factor 
The  residual  gas  in  an  evacuated  ring  will  cause  scattering  losses 
of  the  traversing  beam. 

The  total  scattered  power  Psc  is  calculated,  and  the  quality  factor 
Qsc  is  then  obtained  by 

(W  -  energy  stored)  nj 

^sc  "  TP  =  scattered  power) 
sc 

The  effect  on  the  overall  quality  factor  Q  is  then  given  by  the 
contribution  of  Qsc  to 


I.  j  i_ 
o  r  q,. 


(2) 


Eq.  2  follows  from  eq.  1  under  simplifying  conditions. 
Relation  of  scattered  power  to  Qsc: 


« -  L  - > 


We  will  assume  that  the  scattered  power  per  interval  ax,  aPsc,  is 
proportional  to  the  incident  power  P: 


n  7 


We  will  also  assume  that  the  change  of  P  due  to  scattering  is  negligible 
throughout  the  length  L,  i.e. 

Pcr  «  —  PL  .  (4) 

sc  Ksc 

Here  ksc  =  extinction  length  due  to  scattering. 

The  power  P  in  the  cavity  is  simply  related  to  the  energy  W  stored  in 
the  cavity  through  the  following  considerations: 

Assume  a  monochromatic  plane  field  with  rms  magnitudes  £^^5  and 
Hrm$  in  vacuum.  With  the  Poynting  vector  S  =  E  x  H  =  Uo/uo)  '  .  the 

time-averaged  power  is 

P  =  <  /  §  .  3A  >  =  (e0/u0)1/2  E2rm$A  .  (f-Lff-LX)  (5) 

On  the  other  hand,  the  energy  W  is  given  by 

W»-j/(£.U  +  ff.  ff)dVol  =  e0E2RMSVol  =  e0E2RMSLA 


(cylindrical  mode  volume;  this  assumption  is  however  not  critical). 

We  have  therefore  from  eq.'s  5)  and  6): 

P  =  W/(L/eoITd)  =  (c/l)W  =  (free  spectral  range)x(energy  stored)  (7) 
From  eq's  1),  4),  and  7),  we  get  then 


>sc  =  “0  r 


[L/c]P_ 

(PL/<e. 


We  have  thus  to  increase  the  scattering  length  <sc  such  that  the 
resulting  Qsc  >>  Q.  The  fact  that  we  see  the  sun  indicates  that  this  is 


easy. 


/2<? 


In  the  following,  the  extinction  length  due  to  Rayleigh  scattering 
is  calculated. 

Rayleigh  Scattering 

We  assume  dipole  scattering  only,  because  it  usually  has  the 
largest  scattering  cross  section. 

We  assume  Rayleigh  scattering  off  individual  atoms,  for  two 
reasons: 

1)  The  atomic  dimensions  d  are  small  compared  to  the  wavelength  \  used 
(no  complications  from  charge  distribution  in  atom):  dAIR  =  3.03  x 
10"^  m,  \  =  5  x  10“^  m 

2)  the  mean  free  path  length  i  is  large  compared  to  the  wavelength  used 
(no  scattering  from  density  fluctions) 

l  =  kT/(  /21  7t  d^p)  >>  X  requires  for  air:  p  <<  18  kPa  (=  135  torr) 

We  have  then  for  the  total  scattered  power  off  one  (induced)  dipole 

psc,l  =  (47t3/3Hc/e  MWoA4)  (9) 

where  the  dipole  moment  WQ  =  aE0  =  (9a)  in  an  electric  field  E0, 

with  the  pol arizeabil ity 

ot  =  Uq/NJU  -  1)  -  U0/Nv)(n2  -  1)  (10) 

r\c  =  e  (macroscopic)  relative  dielectric  constant  of  gas, 

Nv  =  volume  density  of  scatterers,  in  an  ideal  gas:  p  =  Nv  kT  (10a). 

(p  =  pressure) 

The  (total)  scattering  cross  section  per  molecule  is 

«  3  2  2  4 

a  s  t  it  (e-l)/Nx, 

O  V 


(11) 


and  the  relation  for  the  total  scattered  power  on  traversing  a  length  L 

is,  using  equations  4,  5,  8,  9,  9a,  10,  and  c/u^- =  1: 

2  2 

=  Per  ,N  AL  =  |  n3  P  =  — PL  ,  (12) 

sc  sc,  1  v  3  ,  < 

’  Nx  sc 

V 

from  which  the  extinction  length  <sc  can  be  calculated,  with  eq.  10a: 

,-1  8  3  (n2-  l)2  32  3  (n  -  1 ) 2kT 

Ksc  '  I  *  =  “1  *  -  . 

N  X  X  p 

v  r 

using  (n^-l)^  =  (n-1 )^(n+l )2*(n-l )^4. 


Using  the  Lorentz-Lorenz  formula  and  eq.  10a: 
2 

n-1  P 

— 2 —  -  n-1  =  const  N  =  const-jTy  ,  or 

n  +1 


(n-1)  =  ( P/P o ) ( T 0/T ) ( n- 1 )  T  .  . 

p.l  .x  u  Pq  » '  OtA 

At  standard  conditions  (p0  =  101  kPa  =  1  atm,  T0  =  273  K  =  0°C),  and  X 
514.5  nm,  nAIR  -  1  =  2.79  x  10'4,  nHe  -  1  =  0.3500  x  10'4. 

We  have  finally 

-1  -  32  3  i  1  T°  p_  . 

Sc  'i'  7*  r  ( ""^po.Tn.x  *  (14) 

where  N0  *  p0/k T0  =  2.687  x  10^5  molecules/m^  (=  Loschmidt  number) 

Quality  factor  versus  scattering: 

The  final  result  is  then,  using  eq.'s  14  and  8): 


Qsc  -  2* 


(sc  .  3  "°  xJ  T  Po 

x  T5  7^r  T7  p 

*  n  T  i 

Po» 1 0»* 


Given  the  type  of  gas  (which  determines  n): 


Example:  1)  Air-filling  at  standard  conditions,  A  =  0.5145  x  10  m 


3 

16 


2.6/8x10  x  (0.5145x10  ) 

— 7 - Tin - 

*  (2.940x10  ) 


8.0x10 


1 1 


Compare  this  to  a  square  ring  with  overall  finesse  F  =  20,000. 

The  latter  results  in  a 

QMirror  =  FL/A  =  1.2  x  1012  (L  =  31  m,  A  =  0.5145  x  10‘6  m) 

which  is  about  equal  to  the  quality  factor  due  to  scattering  by  air  at  1 

atmosphere. 

With  very  modest  forepumping  down  to  1  torr  =  133  Pa,  Qsc  =  6  x  10^ 
i.e.  negligibly  large  compared  to  even  the  best  hoped-for  (achieveable?) 
mirror  quality  factors. 

Note  also  that  the  optimum  pressure  for  HeNe-plasma  for  tubes  of  “  1  cm 
diameter  is  somewhat  below  1  torr;  this  means  that  scattering  losses  are 
negligible  even  in  active  ring  lasers,  if  the  assumptions  above  are 
satisfied  (however,  scattering  in  an  ionized  gas,  and  at  resonating 
atoms  would  have  to  be  considered  separately). 

Example  2:  A  cavity  of  length  L=52  cm  produced  by  Rockwell  has  an 
evacuated  finesse  of  F  =  29800  ±  400  at  A  =  633  nm.  What  is  the 

VaC 

finesse  Fsc  with  air  filling? 

<n-1>A1r,  T„,P0,X  •  633nm=  0-°0029170'  T=(273  -  25)K,  P-P,. 

we  get  Qsc  =  1.661 xlO 1  ,  or  Fsc  =  Qsc A/L  =  2.02xl06. 

From  eq.  2  follows  =  -r-^—  +  +  - - — g- 

Air  vac  rsc  ^  uu  2.02x10 

or  FA.r=  29370 

(measured  at  Rockwell  with  air  filling:  2900  ±  400) 


The  datd  are  compatible  with  our  calculations,  although  they  are  not 


conclusive.  The  same  mirrors,  put  into  a  30  m  cavity,  would  have  a 
finesse  in  vacuum  of  29800x30/0.52  =  1.719x10^,  and  letting  air  in  would 
reduce  the  finesse  from 


6  6 

F  -  1.72x10  to  F..  =  0.929x10  .  ,  which  would  then  mean  a 

V  aC  n I i 

substantial  reduction  in  overfall  finesse,  or  quality  factor. 


III.  Q.  Proposed  calibration  methods  for  large  ring:  Light  drag 

through  gas  flow  from  switchable  calibrated  leaks,  tilting  (rocking) 
the  ring  base. 

The  dynamic  range  of  a  large  stationary  ring  is,  expressed  in  earth 
rotation  rates  from  about  to  or,  expressed  in  optical 

beat  frequencies,  from  about  100  nHz  to  1  kHz,  i.e.  over  10  orders  of 

magnitude.  It  turns  out  that  low  beat  frequencies  can  easily  be 

generated  via  air  flow  through  calibrated  leaks:  The  needed  mass  flow 
requires  such  low  air  pressures  that  the  cavity  quality  factor  (section 
III.  C.)  is  not  affected. 

This  contribution  is  already  published  (ref.'s  20  and  21).  Rocking 

of  the  base  is  made  possible  at  the  base  at  the  Seiler  lab.  It  is 

useful  for  generating  large  signals. 
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Abstract 

Two  types  or*  e  f  f ec  t  s  are  i.nvestigatec  as  :  a  their  potential  '0:  calibrating  3  large 
ring  ana  for  introducing  test  signals  a'  <f,o»n  magnitude:  1)  tilting  the  base  of  the 
ring  oy  a  <no«n  angle  ana  with  <-o*n  t.me  evolution,  ana  2)  utilizing  Fresnel  Crag  'ram 
a  cantroilea  flaw  of  gas. 


Introduction 


A  large,  58  meter  square,  ring  laser  gyro  .RIG)  is  being  assembled  on  the  pneumatic 
isolation  platform  at  tne  ^ran.K  c.  Seiler  Research  .aooratory.  T«o  techniques  far 
calibration  of  the  RLG  will  pe  consioered  in  this  presentation.  In  the##f^rst  method, 
t.ne  RLG  would  oe  tilted  at  a  precisely  known  angle  to  vary  the  ef fective”’fotationai 
rate.  The  second  method  oiscussed  makes  use  of  the  Fresnel  drag  from  a  controlled  flow 
of  gas  in  the  ring. 

An  adjustable  effect  to  produce  an  4lG  output  is  described  for  several  reasons:  1) 
to  check  the  linearity  of  the  output,  2)  to  check  properties  in  the  frequency  domain,  3 
to  calibrate  the  sensitivity  and/or  the  noise  sources  acting  in  the  ring,  and  4)  to 
introduce  a  probe  signal  during  initial  adjustments.  The  earth  rotation  naturally 
provides  a  sizeaole  rotation  rate  with  a  very  good  signai/noise  ratio  of  the  order  of 
1C)3;1  but  it  has  tne  disadvantages  of  not  oeing  adjustable  and  of  being  constant  as 
opposed  to  an  effect  leading  to  a  frequency  modulation  or  the  output.  One  advantage  of 
basing  the  large  RIG  on  a  pneumatically  supported  platform  is  that  the  oase  can  oe 
tilted  by  small  angles  thus  changing  the  effective  rotational  rate  imposed  on  the  ring. 
The  isolation  platform  can  oe  maintained  at  a  known  tilt  to  local  level  to  and  accuracy 
of  the  order  of  10-’  arc  seconds  by  means  of  closed  loop  control  of  pneumatic 
actuators2.  The  platform  can  be  offset  from  level  oy  as  much  as  6  arc  minutes  by 
means  of  these  pneumatic  actuators.  The  base  can  also  be  excited  in  6  degrees  of 
freedom  at  small  amplitudes  from  DC  to  10  Hz  by  means  of  electro-magnetic  actuators. 

The  light  path  of  the  RLG  is  contained  in  tubing  evacuated  to  about  ID*8  torr.  The 
Fresnel  drag  calibration  technique  would  be  performed  by  allowing  a  small  gas  flow 
parallel  to  the  light  beam,  along  a  portion  of  the  ring.  A  Known  gas  flow 

rate  establishes  a  calculable  effect  on  the  difference  in  velocity  for  the  clockwise  an 
counterclockwise  light  beams  of  the  RLG,  and  an  expected  value  for  the  non-rotating 
output  under  this  condition. 


Range  of  calibration 

Large  ring  lasers  are  attractive  because  of  their  high  sensitivity  for  rotation  and 
low  noise.  In  the  following  we  will  use  a  square  ring  laser  of  A  =  58  m2  with 
perimeter  L  *  4  /  58  m,  operating  at  a  wavelength  X  =  514.5  nm  as  basis  for  discussion. 
A  two  mode  ring  laser  has  a  beat  frequency  output 

Af  »  4  a’-Q/A  L,  ( 

where  A  »  surface  vector  of  ring,  “  *  angular  velocity  vector  of  absolute  rinq 
rotation. 
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As  Equation  (1)  is  linear  in  a  ,  there  is  no  obvious  upper  limit  of  output.  However, 
for  convenience  we  will  limit  the  design  considerations  to  Dandwidths  up  to  earth  rate 
for  tne  recording  of  the  phase  t  (t),  [  Af  =  (l/2»)(d<  /dt)],  i.e.,af  =  679  Hz 
(latitude  0  =  33°  58').  As  to  the  lower  limit  we  will  assume  that  quantum  noise 
limits  the  performance  at  high  Fourier  frequencies.  The  power  spectral  density  is  given 
Dy  8 


Sdf,w  •  2  <^o/Q2?) 

with  an  estimate  of  finesse  F  =  1000  or  a  quality  factor  Q  =  (L/l  )F  =  5.9  x  10^,  ana 
a  power  level  of  500  mw,  the  white  power  spectral  density  would  be,  at  the  oscillation 
frequency  f0  =  c/X  =  583  THz , 

S  4f,»  =  1.51  X  10-10  (Hz2/H2) , 

which  allows  an  accuracy  of 


5Af  =  ^SOf,w/ 


(3) 


for  the  measurement  of  the  beat  frequency  over  an  integration  time  r;  e.g.,  with  r  = 

1  o,  an  rms  frequency  fluctuation  of  the  order  of  42  nHz  can  oe  anticipated.  However, 
the  open-enoed  limit  set  oy  Equation  (3)  is  not  realizable  over  arbitrarily  long  times. 
In  practice,  1/f  noise  or  even  l/f2  noise  of  the  laser  oscillators  set  a  lower 
limit4  to  the  attainable  accuracy.  In  today's  commercial  ring  laser  gyros,  the  white 
quantum  noise  turns  over  into  1/f  noise  at  corner  frequencies  of  the  order  of  10-1000 k 
hz.  Furthermore,  tne  observed  1/f  noise  is  always  larger  than  the  Gagnepain-Uebersfeld 
limit  as  observed  in  quartz  oscillators^. 


SAf.F  *  A/f>  A  *  fo/(34  (A) 

For  lack  of  oetter  theory,  we  use  Equation  (4)  to  estimate  the  flicker  floor,  with  f0 
=  583  THz  and  Q  =  5.9  x  10l°,  S  if,F  *  2.8  x  10"i4/f.  The  equivalent  Allan 
standard  deviation 

o,(Tj  =  v'2in2  x  2.0X1O-14  =  197  nHz. 

The  transition  between  white  noise  and  1/f  noise  appears  at  Fourier  frequencies  of  about 
200  i*Hz,  requiring  integration  times  of  about  1  hour  to  arrive  at  this  flicker  floor; 
lower  limits  of  oeat  frequency  errors  of  the  order  of  10~8  Hz  are  theoretically 
feasible,  and  a  dynamic  range  of  over  nine  decades  of  beat  frequency  is  possible,  for 
wnicn  range  calibration  signals  should  be  provided. 

First  calibration  method:  tilting  the  gyro  base 

A  tilt  modifies  the  output  from  earth  rate.  The  surface^T  in  Equation  (1)  is  made 
time-dependent,  A  — A(t).  This  can  be  achieved  by  rotating  an  actively  controlled  pad 
around  the  east-west  symmetry  axis.  A  variety  of  factors  have  to  be  investigated  in 
such  a  feasibility  study;  Reproducibility,  adjustibility ,  hysteresis,  effect  on 
long-term  stability.  The  change  of  beat  frequency  due  to  tilt  in  the  meridional  plane 
by  an  angle  ia  is 


dAf(Hz)  *  (4  A  •  5/ A  L  )  cot  9  6a  *4. lx  10“^  8a  (arc  second)  ( 

in  the  example  above.  To  change  the  beat  frequency  by  1  mHz,  a  lift  of  Qj  u m  of  the 
7.62  m  -  long  pad  is  required.  The  "noise"  of  the  tilt  is  10-3  arc  second  which  is 
transduced  by  the  ring  into  a  beat  frequency  fluctuation  of  6.58  uHz,  well  above  the 
limit  given  oy  the  flicker  floor. 


There  is  also  the  possibility  of  rocking  tne  pad,  i.e.,  to  introduce  a  sinusoidal 
.lit  versus  time.  Provided  that  the  excitation  of  eigenmodes  of  the  pad  can  be 
sufficiently  suppressed,  this  manner  of  introduction  of  an  ac-modulated  frequency 
becomes  an  excellent  means  of  generating  a  calibration  signal.  An  external  measurement 
or  uhe  neignt  by  using  high-performance  levelling  instruments  Oy  electrical  or  optical 
means  can  certainly  oe  done  with  an  error  of  10  yw>  and  possibly  of  1  pm. 
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Second  calibration  metnoa 


Light  drag  is  a  surprisingly  simple  means  to  introduce  an  anisotropy  in  the  ring. 
Experimentally,  some  drag  equations  nave  oeen  verified  to  me  level  of  3. IX.6 

for  a  homogeneous  velocity  fLeid  of  a  gas  parallel  to  the  light  oeam  over  a 
distance  d,  the  beat  frequency  is  to  a  good  approximation  given  by7 

Af  s  (4/ X  L)  (n-1)  vd,  (6 

„nere  n  =  index  of  refraction.  For  helium  this  npproximation  is  good  to  IX.  The 
quantity  (n-1)  is  proportional  to  the  gas  pressure.  For  an  ideal  gas 


”  (VT)(p/Po)(n-l)TojPoiX  '7 

T  =  aosolute  temperature  of  flowing  gas,  p  =  pressure  of  flowing  gas,  TJ(  p0  = 
standard  conditions.  For  stationary  mass  flow,  the  flow  rate 

m  =  d(pv)/dt  =  q(vp)  *  kT  dN/dt  »  constant  (8 

(q  =  tube  cross  section,  V  =  volume  of  gas,  N  =  number  of  gas  molecules,  k  =  Boltzman 
constant),  if  a  mass  flow  source  is  used,  for  example  a  calibrated  leak  as  is  used  in 
vacuum  work.  It  follows  for  the  beat  frequency  due  to  light  drag 


Af  -  (4/XL) (d/q) (n-l)T  -  (T0/Tp0)M  '9 

O  t rO  f A 

The  beat  frequency  is  therefore  determined  by  the  mass  flow  rate.  The  velocity  vs 
iength  is  arbitrary  and  need  not  be  known;  the  drag  is  Independent  of  pressure  and 
pressure  gradients  along  the  path,  as  long  as  M  is  constant  and  the  flow  is  stationary. 

Estimate  of  required  leak.  Assume  a  helium-filled  leak  with  (n-l)T  p_^  =  3.49 
x  1CT3.  The  leak  is  to  be  injected  into  one  leg  of  the  ring  having°l  Cm  tube  diameter 
and  d  = ,J.?J.ength .  At  room  temperature  (T  *  298  K),  a  flow  rate  of  M  «2,C  x  10-4 
Pa.m3/s  is  then  required  to  produce  4  f  *  1  mHz.  A  calibrated  He  leak,  of  100  "micron 
cubic  foot/hour"  has  a  mass  flow  rate  of  1.05  x  10"4  Pa.m^/s.  If  larger  beat 
frequencies  are  required,  air  can  be  used  with  about  eightfold  increased  drag. 

Several  calibrated  leaks  can  be  mounted  in  opposing  legs  of  the  ring,  to  change 
magnitude  and  sign  of  the  drag.  They  can  be  switched  off  and  on  at  any  time. 
Commercially,  a  large  range  of  10-11  Pa«-m3/s  to  10~4  Pa«m3/s  is  available  which 
will  cover  the  entire  range  of  beat  frequencies  from  about  10  mHz  down  into  the  noise. 

It  may  also  be  mentioned  that  turbulent  flow  will  remove  the  radial  velocity  profile  and 
will  produce  a  more  uniform  output  frequency.8 

Vacuum  acquirements.  Conversely,  inadvertent  drag  due  to  moving  residual  gas  has  to  be 
avoided. From  the  foregoing  we  can  estimate  upper  limits  of  mass  flow  sources  for 
negligible  drag.  To  push  drag  effect  below  10“7  Hz  generally  we  require  the  sources 
(inadvertent  leaks,  outgassing,  etc.)  to  have  flow  rates  below  10-9  Pa»m3/s. 

However,  typical  information  about  vacuum  is  given  through  pressure,  which  in  turn  is 
related  to  flow  via  the  pump  resistance  of  the  ring. 

It  turns  out  that  with  mass  flow  rates  even  as  high  as  100  micron  cubic  foot  per  hour 
(10'*  Pa«md/s)  we  are  already  in  the  high  vacuum  regime.  In  this  regime,  the  pump 
resistance  R  (s/m3)  of  a  cylindrical  tube  of  cross  section  q  and  length  d,  molecular 
weight  m  is9 


R(s/m3)  -  (3»/4) (m/2kT) 1/2  d/q3/2 

and  the  pressure  difference  along  such  a  tube  is 


(10) 


P2  -  P1  -  HR. 


(11) 
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V  , 


Assuming  a  worst  case  of  "inacvertent"  mass  flow  injection  of  N  =  10-9  Pa  iJ/s 
farthest  away  from  tne  pumo ,  a  =  7  .(■  m,  ana  a  tube  with  1  cm  diameter  as  above,  we  ootain 

P2  .  P,  rp2  =  Pa  ( -h  7  x  10-7  torr) 

wnereoy  tne  pressure  Pi  at  the  input  of  a  75  1/s  pump  would  be  1.3xi0'8  Pa 
1*10-10  torr),  i.e.,  negligible. 

Prom  this  we  conclude  that  reasonable  nigh  vacuum  conditions  have  to  prevail,  but 
also  that  the  loading  of  vacuum  pumps  by  Fresnel  drag  calibrations  is  not  excessive. 
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